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Optical  Properties  of  the  Atmosphere 

(Third  Edition) 

1.  IMRIIIM  .'.'HUN 

Tins  report,  is  intended  to  provide  data  from  which  the  attenuation  or  emission 
of  radiation  along  a  path  in  a  given  model  atmosphere  can  be  determined,  The 
sptcuol  region  covt  n  i  is  from  0.25  to  '.15  micrometers.  We  have  not  concerned 
ourselves  with  the  atmosphere  above  100  kin,  thereby  neglecting  less  than  one- 
millionth  of  the  total  atmospheric  mass. 

We  have  defined  five  model  atmospheres  each  for  temperature,  pressure  and 
absorbing  gas  concentrations.  In  addition  to  these,  two  aerosol  models  have  been 
defined.  Thus,  one  effectively  has  10  models  from  which  to  choose  in  order  to  use 
the  data  as  presented  without  interpolation.  However,  interpolation  senemes  are 
descr  ibed  in  or  der  to  allow  the  user  to  obtain  optical  data  for  intermediate  condi¬ 
tions. 

As  it  is  not  possible  to  cover  all  conceivable  atmospheric  models,  neither  is 
it  possible  to  provide  data  in  a  short  report  for  any  spectral  resolution.  Thus,  we 
are  forced  to  limit  !he  high  resolution  calculations  to  a  few  specific  wavelengths. 

A  list  of  laser  emission  wavelengths  was  chosen  for  this  purpose.  We  attempted  to 
pick  those  lasers  most  commonly  used  in  practical  applications  and  as  research 
tools.  In  addition  to  the  laser  beam  transmittance  table,  we  have  provided  in 
graphical  form  a  calculation  scheme  for  determining  atmospheric  absorption  and 

(Received  for  publication  22  August  1072) 
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scattering  at  about  a  20  wavenumber  resolution  across  the  entire  spectral  region. 
These  curves,  although  appropriate  foi  low  resolution  transmittance  calculations 
would  produce  serious  error  if  used  for  estimates  of  laser  propagation  in  the 
atmosphere.  For  high  resolution  curves  appropriate  to  the  laser  propagation 
problem,  additional  reports  have  been  published  (McClaichey,  197]  and  McClatchey 
&  Selby,  1072). 

The  models  defined  in  Section  2  together  with  the  low  resolution  absorption 
calculation  techniques  presented  in  Section  5,  were  used  to  compute  the  radiation 
emitted  to  space  and  the  radiation  as  observed  looking  up  at  the  sky  from  the  sur¬ 
face  (sky  background)  for  wavelengths  longer  than  5  jmi  (Section  (i).  The  radiation 
to  space  and  sky  backgrounds  at  shorter  wavelengths  is  a  complicated  function  of 
the  solar  zenith  angle  and  the  angle  of  observation,  so  only  a  few  specific  examples 
have  been  given. 

The  effect  of  the  refractive  index  on  the  geometrical  path  length  and  the  effect 
of  cloud  and  fog  on  the  attenuation  of  radiation  are  discussed  in  Sections  7  and  B, 
respectively. 

2.  MODIll.  UMOSIMlkHUS  I  SKI)  IN  <11*  1 U.  U.  Cil.l  l  l.VHONS 

2.1  Standard  Madid  Vtiiiojqihrrej- 

A  series  of  five  model  atmospheres  is  presented  in  Table  1.  The  pressure, 
temperature,  density,  water  vapor  concentrations  are  those  of  the  U.S.  Standard 
Atmosphere  of  1932  and  the  Supplement  Atmospheres  as  compiled  in  the  Hand¬ 
book  of  Geophysics  and  Space  Environment  (Valley  1965).  The  open.*  concentra¬ 
tions  were  extracted  from  Chapter  6  of  the  Handbook  of  Geophysics  and  Space  En¬ 
vironment  (ValLcy,  1969).  The  water  vapor  densities  above  11  km  have  been  taken 
from  Sissenwine  et  al  (1908).  A  discussion  of  units  and  conversion  can  be  found 
in  Appendix  A. 

2.2  t'lonrcnttutions  iif  hinfonnh  Mixed 

Table  2  is  intended  to  provide  the  reader  with  the  data  for  determining  the 
concentrations  of  all  uniformly  mixed  atmospheric  gases  in  any  horizontal  or  ver¬ 
tical  path.  Column  3  gives  the  concentrations  in  parts  per  million  by  volume.  The 
fourth  column  gives  the  total  amount  of  gas  in  a  vertical  path  through  the  entire 
atmosphere.  The  fifth  column  lists  the  numoer  of  (cm -atm of  gas  in  a  one 
kilometer  path  at  sea  level.  The  number  of  (cm  -atm)g^.p  in  a  horizontal  path  at 
any  other  altitude  is  given  by  Eq.  (1). 
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Table  1.  Model  Atmospheres  Used  as  a  Basis  of  the  Computation 
of  Atmospheric  Optical  Properties 


TROPIC'  A  L. 

lit. 

Pressure 

Temp. 

Density 

Water  Vapor 

Ozone 

(km) 

(nib) 

(CK) 

(g/rr.3) 

(g/m3) 

(g/m3) 

0 

1.  013E+03 

300.  0 

1.  1G7E+03 

1. 9E+01 

5.  GP' -05 

1 

9.  04  OK +  02 

294.  0 

1.  OG4E+03 

1.  3E+01 

5.  GE-05 

2 

8.  050K+02 

288.  0 

9. G89E+02 

9.  3E+00 

5.  4E-05 

3 

7.  13  OK +02 

284.  0 

8. 75GE+02 

4. 7E+00 

5.  IE -05 

4 

0.  330K+02 

277.  0 

7  951E+02 

2.2E+00 

4.7E-05 

J 

5. 590E+02 

270.  0 

7.  199E+02 

1.  5E+00 

4.  5E-05 

G 

4. 920E+02 

2G4.  0 

6. 501E+02 

8.  5E-01 

4. 3E-05 

4.  32 0E +02 

257. 0 

5. CS5E+02 

4.  7E-01 

4.  1E-C5 

H 

3.780K+02 

250.  0 

5.  2  L)8E  +  U2 

2.  SE-01 

3. 9E-05 

9 

3. 290K+02 

244. 0 

4. 708E+02 

1.  2K-01 

3.9E-0S 

10 

2. 860E+02 

237. 0 

4.  202E+02 

5.  0E -02 

3.  9E-05 

11 

2.  47  0E +  02 

230.  0 

3.  740E+02 

1.  7K-02 

4.  IE -05 

12 

2. 130E+02 

224. 0 

3.  316E+02 

G.  0E -03 

4.  3K-05 

13 

1.  02  OK +02 

2  17.0 

2.  929K+02 

i.  8PJ  -  03 

4.  GE-05 

14 

t.  5G0E+02 

210.  0 

2. 578E+02 

1.  0E -03 

4  5E-05 

15 

1.  32  OK +02 

204.  0 

2. 20 op: +02 

7.  GK-04 

4.  7E-05 

10 

1.  110K+02 

197.  0 

1.  972E+02 

6.  4K-04 

4.  7E-05 

17 

9.  370E+01 

195.  0 

1.  G7GK+02 

5.  GP: -04 

0.  9E-05 

18 

7. 890E+01 

169.  0 

1.  3B2E:+02 

5.  0E  -04 

9-  0E-05 

19 

fi. 6G0E+01 

203.  0 

1.  14  5P>02 

4.  9E  04 

1.  4E-04 

20 

5. 650E+01 

207.  0 

9.  515E+01 

4.  5E-04 

1.  9E-04 

21 

4. 800E+01 

211.  0 

7 . 938E+01 

5.  IE -04 

2.  4E-04 

22 

4. 090E+01 

2  13.  0 

G. 649E+01 

5.  IE -04 

2.  8E-04 

23 

3. 500E+01 

2  17.  0 

5. GieE+01 

5.  4E-04 

3-  ZE-G4 

24 

3. 000E+01 

2  19.  0 

4. 7G3E+01 

G.  0E-04 

3.  4E-04 

25 

2. 570E+01 

221.  0 

4. 045E+01 

G.  7E-04 

3.  4E-04 

30 

1.22  0E+01 

232.  0 

1.  831E+01 

3.  6E-04 

2.  4E-04 

35 

6.  OOOE+OO 

243.  0 

8.  600E+00 

1.  IE  -04 

9.  2E-05 

40 

3.  050E+00 

254.  9 

4.  181E+00 

4.  3E-05 

4.  IE-05 

45 

1.  590E+00 

2G5.  0 

2.  097t:+00 

1.  9E-05 

1.  3E-05 

50 

8.  540E-01 

270.  0 

1.  101E+00 

6.  3E-06 

4.  3E-06 

70 

5. 790E-02 

219.9 

9.  2  10E-02 

1.4E-07 

8.  6E-08 

100 

3.  000E-04 

2  10.  0 

5.  000E-O4 

1.0E-00 

4.  3E-11 

* 

I 
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Table  1  (Contd).  Model  Atmospheres  Used  as  a  lias  is  of  the  Computation 
of  Atmospheric  Optical  Properties 


MIDI  ATI  TUDli  SUM  MEIf 


lit 

(km) 

Pressure 

(mb) 

Temp. 

(°K) 

Density 

(g/mT) 

Water  Vapor 
(e/m3) 

Ozone 
(g/  m  3 ) 

0 

1.  013E+03 

294. 0 

1.  19110+03 

1. 4K+01 

6.  OE -05 

1 

9.  02010+02 

2  90.  0 

1.  0801  +03 

9.  310+00 

G.  0E -05 

2 

8.  020E+02 

285.  0 

9.7  5710  +  02 

5.  OE-’-OO 

6.  OE -05 

3 

7.  1O0E+02 

27  9.  0 

8.  84610+02 

3.  310+00 

6.  210-05 

4 

6.  28010+02 

273.  0 

7. 998E+02 

1.  910+00 

6.  410-05 

5 

5.  54  010+02 

267. 0 

7.2  1  110+02 

1.  OE+00 

6.  GE-05 

G 

4.  87  OK+02 

261.  0 

6. 487K+02 

6.  110-01 

G.  910-05 

7 

4. 200K+02 

255.  0 

5. 830E+02 

3.7E-01 

7.  5E-05 

8 

3.  72010+02 

24  8.  0 

5.225E+02 

2.  110-01 

7.  9E-05 

9 

3.  24010+02 

242.  0 

4. 66910+02 

1. 210-01 

8.  6E-05 

10 

2.  81010+02 

2  35.  0 

4.  15  OK +02 

6.  410-02 

9.  0E  -05 

11 

2.  43010+02 

229.  0 

3.  69310+02 

2.210-02 

1.  IE  -  04 

12 

2.  09010  +  02 

222.  0 

3.26910+02 

6.  OE -03 

1.  2  E  -  04 

13 

1.79010!  02 

2  16.  0 

2.  8  82!'' +09 

1.  HE  -  03 

1.  5E-04 

l'i 

1. 5  3  OK +02 

2  16.  0 

2.  46410  +  02 

1.  0E -03 

1.  8E-04 

15 

1.  30010+02 

2  16.  0 

2.  10410+02 

7  .  OK  -  04 

1.  9E-04 

16 

1.  110E+02 

2  16.  0 

1. 79710+02 

6.  4E-04 

2.  IE -04 

17 

9.  50010+01 

2  16.  0 

1.  53510  +  02 

5.  OK -04 

2 . 4 10  -  04 

18 

8.  12010  +  01 

216.  0 

1.  305E+02 

5.  OE -04 

2.  8E-04 

19 

6.  950K+01 

2  17. 0 

1.  1  10F+02 

4.  9E-04 

3.  2E-04 

2  0 

5. 950E+C1 

2  18.  0 

9.  4  5310+01 

4.  5E-04 

3.  4E-04 

21 

5. 100E+01 

2  19.  0 

8.  05GE+01 

5.  IE -04 

3.  6E-04 

22 

4. 370E+01 

220.0 

6.  872E+01 

5.  IE -04 

3.  6E-04 

2  3 

3. 760E+01 

222.0 

c  ^  357j^j-0  i 

5.  4E-04 

3.  4E-04 

24 

3. 220E+01 

223.  0 

5.  014E+01 

6.  OE-04 

3.2E-04 

25 

2. 770E+01 

224.  0 

4.  288K+01 

6.7E-04 

3.  0E-04 

30 

1.  320E+01 

234.0 

1.  32210+01 

3.  CE-04 

2-  0E-04 

35 

6.  520E+00 

245.  0 

6.  5I9E+00 

1.  IE-04 

9.  2E-05 

40 

3. 330E+00 

258.  0 

3.  330E+00 

4.  3E-05 

4.  IE -05 

45 

1.  760E+00 

270.  0 

1.757E+00 

1.  9E-05 

1.  3E-05 

50 

9. 510E-01 

276.  0 

9.  512E-01 

C.  3E-06 

4.  3E-0G 

70 

6. 7 10E-02 

218.  0 

6.  706C-02 

1.4E-07 

8.  6E-08 

100 

3. 000E-U4 

210.  0 

5.  OOOE-04 

1.  0E-09 

4.  3E-11 

Table  1  (Contd).  Model  Atmospheres  Used  as  a  Basis  of  the  Computation 
of  Atmospheric  Optical  Properties 


MIDLATITUDE  WINTER 


Ht. 

(km) 

Pressure 

(mb) 

Temp. 

<°K) 

Density 

(g/m3J 

Water  Vapor 
(g/m3) 

Ozone 

(g/m3) 

0 

1.  018E+03 

272.2 

1.  301E+03 

3.  5E+00 

6.  0E-05 

1 

8.  97  3E+02- 

268.7 

1.  162E+03 

2.  5E+00 

5.  4E-05 

2 

7.  897E+02 

205. 2 

1.  037E+03 

1.  8E+00 

4.  9E-05 

3 

6,  938E+02 

261.7 

9.  230E+02 

1.  2E+00 

4.  9E-05 

4 

6. 081E+02 

255.7 

8.  2  82E  *-92 

6.  6E-01 

4.  9E-05 

5 

5.  313E+02 

249-7 

7.41 1E+02 

3.  8E-0) 

5.  8E-05 

6 

4  627E+02 

243.  7 

6.614E+02 

2.  IE -0 1 

6.  4E-05 

7 

4.  016E+02 

237.7 

5.  886E+02 

8.  5E-02 

7.  7E-05 

a 

3.  47  3E+02 

231.  7 

5. 222E+02 

3.  5E-02 

S.  0E-05 

9 

2.  992E+02 

225.  7 

4.  619E+02 

1.  6E-02 

1. 2E-04 

10 

2.  568E+02 

2  19.7 

4.  072E+02 

7.  5E-03 

1.  6E-04 

11 

2.  199E+02 

2  19.  2 

3. 496E+02 

6.  9E-03 

2.  IF. -04 

12 

1.  882E+02 

218.7 

2.  999E+02 

6.  0E-03 

2.  RE-04 

13 

1.  610E+02 

218.  2 

2.  572E+02 

1.  8E-03 

3.  0E-04 

14 

1.  378E+0'-' 

2  17.  7 

2. 2O6Ei-02 

1.  0E-03 

3.  2F-04 

15 

1.  17  8E+0r 

217.2 

1.  890E+02 

7.  6E-04 

3.4E-04 

10 

1.  007E+02 

2  16.  7 

1.  820E+02 

6.  4E-04 

3.  OE-04 

17 

8. 6 10E+01 

216.2 

1.  388E+02 

5.  6E-04 

3.  9E-04 

18 

7.  350E+01 

215.  7 

1.  188E+02 

5.  0E-04 

4.  IK-04 

19 

6. 280E+01 

215.2 

1.  017E+02 

4.  9E-04 

4.  3E-04 

20 

5.  370E+01 

215.2 

8. 690E+01 

4.  5E-04 

4. 5E-04 

21 

4. 580E+01 

215.2 

7.421E+01 

5.  IE-04 

4.  3E-04 

22 

3.  910E+01 

215.2 

0. 338E+01 

5.  IE-04 

4.  3E-04 

23 

3. 340E+01 

215.2 

5. 415E+01 

5.  4E-04 

3.  9E-04 

24 

2. 860E401 

215.2 

4. 624E+01 

0.  0E-04 

3.  6E-04 

25 

2. 430E+01 

2  15.2 

3.  950E+01 

6.  7E-04 

3.  4E-04 

30 

1.  110E+01 

217.4 

1.  783E+01 

3.  6E-04 

1.  9E-04 

35 

5.  189E+00 

227.  8 

7. 924E+00 

1.  IE -04 

9.2E-05 

40 

2. 530E+0C 

243.2 

3.  62  5E+00 

4.  3E-05 

4.  IE -05 

45 

1.290E+00 

258.  5 

1.74 1E+00 

1.  9E-05 

1.  3E-05 

50 

6  820E-01 

265.  7 

8.954E-01 

6.  3E-06 

4.  3E-06 

70 

4. 670E-02 

230.7 

7.051E-02 

1.4E-07 

8.  6E-O8 

100 

3. 000E-04 

210.  2 

5.  000E-04 

1.  0E-09 

4.3E-11 
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Table  1  (Contd).  Model  Atmospheres  Used  as  a  Basts  ol  the  Computation 
of  Atmospheric  Optical  Properties 


SUBARCTIC  SUMMER 

lit. 

(km) 

Pressure 

(mb) 

Temp. 

<°K) 

Density 

(g/m3) 

Water  Vapor 
(g/m5) 

Ozone 

(g/m3) 

0 

1.  010E+03 

287.  0 

1.  220E+03 

9.  1E+00 

4.  9E-05 

1 

8.  960E+02 

282.  O' 

1.  110E-03 

6.  OE  WO 

5.4E-05 

2 

7. 929E+02 

276.  0 

9.  971E+02 

4.2E+00 

5.  6E-05 

3 

7  O00E+O2 

271.  0 

8.  985E+02 

2.7E+00 

5.  8E-05 

4 

6. 160E+02 

266.  0 

8.  077E+02 

1.  7E+00 

6.  0E-05 

5 

5. 410E+02 

260.  0 

7. 244E+C2 

1.  GE+00 

6.4E-05 

6 

4. 730E+02 

253.  0 

6.  519E+02 

5.  4E-01 

7. IE-05 

1 

4. 130E+02 

246.0 

5.  849E+02 

2.  9E-01 

7. 5E-05 

8 

3. 590E+02 

239-  0 

5. 231E+02 

1.  3E-02 

7. 9E-05 

9 

3.  107E+02 

232.  0 

4. 663E+02 

4.  2E-02 

1.  IE -04 

10 

2.  677E+02 

225.0 

4.  142E+02 

1.  5E-02 

1.  3E-04 

11 

2.  300E+02 

225.  0 

3.  559E 1 02 

9.  4E-03 

1.  8E-04 

13 

1.  977E+02 

225.0 

3.059E+02 

6.  0E-O3 

2.  IE-04 

13 

1.  7  00E1-02 

225.0 

2.  t>30Ei-02 

2.  3E-03 

2,  KE-04 

14 

1.460E+02 

225.  0 

2.260E+02 

1.  OE-03 

2.  8E-04 

15 

1.  250E+02 

225.0 

1.  943E+02 

7.  6E-04 

3.  2E-04 

16 

1. 080E+02 

225.0 

1.67  1E+02 

6.4E-04 

3.  4E-04 

17 

9. 280E+01 

225.  0 

1.436E+02 

5.  6E-04 

3.  9E-04 

18 

7. 980E+01 

225.  0 

1.235E+02 

5.  0E-04 

4.  E-04 

19 

6.  860E+01 

225.  0 

1.  062E+02 

4.  9E-04 

4.  IE -04 

20 

5.  890E+01 

22  5.  0 

9.  128E+01 

4.  5E-04 

3.9E-04 

21 

5.  07  0E+01 

225.  0 

7. 849E+01 

5.  IE -04 

3.  6E-04 

22 

4.  360E+01 

225.  0 

6. 750E+01 

5.  IE -04 

3.2E-04 

23 

3.750E+01 

225.  0 

5.  805E+01 

5. 4E-04 

3.  OE-0'i 

24 

3. 227E+01 

226.  0 

4.  963E+01 

6.  0E-04 

2.  8E-04 

25 

2.  780E+01 

228.  0 

4.  247E+01 

6.  7E-04 

2.  6E-04 

30 

1.  340E+01 

235.0 

1.  338E+01 

3.  6E-C4 

1.  4E-04 

35 

6.  G10E+00 

2t7. 0 

6.  614E+00 

1.  IE-04 

9.2E-05 

40 

3.  400E+00 

262.0 

3.  404E+00 

4.  3E-05 

4.  IE-05 

45 

1.  810E+00 

274.  0 

1.817E+00 

1.  9E-05 

1.  3E-05 

50 

9.  870E-01 

277. 0 

9.  863E-01 

6.  3E-06 

4.  3E-06 

70 

7. 070E-02 

216.  0 

7. 07  IE -02 

1.4E-07 

8.  6E-08 

100 

3.  OOOE-04 

210.0 

f>.  O00E-04 

1.  OE-09 

4.  3E-11 
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Table  1  (Contd) ,  Model  Atmospheres  Used  as  a  Basis  of  the  Computation 
of  Atmospheric  Optical  Properties 


SUBARCTIC  WINTER  | 

Ht. 

(km) 

Pressure 

(mb) 

Temp. 

<°K) 

Density 

(g/mJ) 

Water  Vapor 
(g/m3) 

Ozone 

(g/m3) 

0 

1.  013E+03 

257.  1 

1.  372E+03 

1.  2E+00 

4.  IE-05 

1 

8.  87  8E+02 

259.  1 

1. 193E+03 

1. 2E+00 

4. IE-05 

2 

7. 775E+02 

255.  9 

1.  058E+03 

9.  4E-01 

4.  IE-05 

3 

fi.  798E+02 

252.7 

9,  366E+02 

6.  8E-01 

4.  3E-05 

4 

5.932K+02 

247.7 

8.  339E+02 

4. IE-01 

4.  5E-05 

5 

5.  158F.+02 

240.  9 

7 . 457E+02 

2. 0E-01 

4.7E-05 

6 

4. 467E+02 

234.  1 

6. S46E+02 

9.  8E-02 

4.  9E-05 

7 

3.  85  3E+G2 

22  7.  3 

5.  904E+02 

5.4E-02 

7.  IE-05 

8 

3.  308E+02 

220.  6 

5. 226E+02 

1.  IE -02 

9.  GE-05 

9 

2. 829E+02 

2  17.  2 

4. 538E+02 

8.  4E-03 

1.  6F-04 

10 

2.418E+02 

217.2 

3.  879E+02 

5.  5E-03 

2.  4E-  04 

11 

2.  067E+02 

217.2 

3.  315E+02 

3.  8E-03 

3.2E-04 

12 

1.  766F+02 

217.2 

2. 834E+02 

2.  6E-03 

4.  3E-04 

13 

1.  510E+02 

2  17.2 

2. 422E+02 

1.  Gib  - 03 

4.  7E-04 

14 

1.  291E+02 

217.2 

2.07 1E+02 

1.  0E-03 

4.  9E-04 

15 

1.  103E+02 

2  17 . 2 

1.  770E+02 

7.  6E-04 

5.  6E-04 

16 

9. 431E+01 

216.6 

1.  517E+02 

6.  4E-04 

6.2E-04 

17 

8. 058E+01 

216.0 

1.  300E+02 

5.  6E-0^ 

6.2E-04 

18 

6.  882E+01 

215.4 

1.  1 13E+02 

5.  0E-04 

6.2E-04 

19 

5.  875E+01 

214.8 

9.  529E+01 

4.  9E-04 

6.  0E-04 

20 

5.  014E+01 

214.  1 

8.  155E+01 

4.  5E-04 

5.  6E-04 

21 

4. 277E+01 

213.6 

6. 976E+01 

5.  IE -04 

5.  IE -04 

22 

3.  647E+01 

213.0 

5. 966E+01 

5.  IE-04 

4.  7E-04 

23 

3. 109E+01 

2  12.4 

5.  100E-*01 

5.  4E-04 

4.  3E  -04 

24 

2. 649E+01 

211.  8 

4.  358E+01 

6.  0E  -04 

3.  6E-04 

25 

2.256EHJ1 

2  11.  2 

3. 722E+01 

6.  7F-04 

3.  2E-04 

30 

1.  020E+01 

216.0 

1. 645E+01 

3.6E-04 

1.  5E-04 

35 

4.701E+00 

222.  2 

7. 368E+00 

1.  IE-04 

9.2E-05 

40 

2.243E+00 

234.7 

3. 330E+00 

4.3E-05 

4.  IE-05 

45 

1.  113E+00 

247.  0 

1.  569E+00 

1.9E-05 

1.3E-05 

50 

5.719E-01 

259.  3 

7. 682E-01 

6-  3E-06 

4.3E-06 

70 

4. 016E-02 

245.7 

5.  695E-02 

1.4E-07 

8.  6E-08 

100 

3.  00QE-04 

2  10.  0 

5.  000E-04 

1.  0E-09 

4.3E-11 

Table  2.  Concentrations  of  Uniformly  Mixed  Gases 


A 


Constitu¬ 

ent 

Mole¬ 

cular 

V/t. 

ppm  by 
Vol. 

(cm-atm)g-pp 

in  vertical  path 
from  sea  level 

(cm-atrn)STP/km 
in  horizontal 
path  at  sea  level 

gm  cm  2/mb 

lief. 

Air 

26.  97 

106 

8xl05 

105 

1.  02 

* 

co2 

44 

330 

264 

33 

5.  llxlO'4 

t 

N2o 

44 

0.  28 

0.22 

0-  028 

4.  34xl0'7 

+ 

CO 

28 

0.  075 

O.OS 

0.0075 

7. 39x10'® 

0 

ch4 

16 

1.  6 

1.28 

0-  16 

9.  OlxlO"7 

** 

°2 

32 

2. 095xl05 

1.  68xl05 

2. 095xl04 

0.  236 

* 

where  L  is  the  path  length  in  km,  and  PQ  and  Tq  refer  to  STP  (PQ  =  1  atm  and 

T  =  273K). 
o 

Tile  sixth  column  provides  values  from  which  the  amount  ol  gas  between  any 
two  pressure  levels  in  the  atmosphere  can  be  determined  by  multiplication  of  the 
tabulated  values  by  the  pressure  increment  in  millibars  (Ap)  times  the  secant  of 
the  zenith  angle  (see  Appendix). 

2.3  Aerosol  Oislti  billions 

The  two  aerosol  models  (Table  3)  describe  a  "clear"  and  a  "hazy"  atmosphere 
corresponding  to  a  visibility  of  23  and  5  km  respectively  at  ground  level.  The 
aerosol  size  distribution  function  for  both  models  is  the  same  at  all  altitudes  and 
is  similar  to  the  one  suggested  by  Delrmendjian  (1963)  for  continental  haze 
(Figure  1).  whereby  the  relationship  between  the  number  density  [n(r)  J  and  the 
particle  radius  is  given  by: 

— --  =  C.  -  r  "*  for  0.  1  pm  <  r  <  10  pm 
Ar  1 

«  Cj*  104  for  0.  02  (im  <  r  <  3.  1  pm  (2) 

*  =0  for  r  <  0.  02  and  r  >  10.  0. 

Ar 

It  differs  from  Deirniendjian's  model  "C"  in  that  the  large  particle  cutoff  has  been 
extended  from  5  pm  to  10  pm.  This  distribution  function  is  normalized  with 

•Valley  (1965),  f  Fink  et  al  (1964)  ,*Birklar.d  and  Shaw  (1959),  °Shaw  (1968), 
**5haw  (1969). 
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Table  3.  Aerosol  Models  -  Vertical  Distributions  for  a  "Clear" 
and  "Hazy"  Atmosphere 


PARTICLE  DENSITY  „ 

N  (PARTICLES  PER  crn  ) 

Altitude 

23  km  Visibility 

5  km  Visibility 

(km) 

Clear 

Hazy 

0 

2. 828E+03 

1.  378E+04 

1 

1.  244E+03 

5. 030E+03 

2 

5.  37  1E+02 

1 . 844E+03 

3 

2.  25RE+02 

6.731E+02 

4 

1.  192E+02 

2. 453E+02 

5 

8.  987E+01 

8.  987E+01 

6 

G.  337E+01 

6.  337E+01 

7 

5.  890E+01 

5. 890E+01 

8 

6.  069E+01 

6.  069E+01 

9 

5.  818E+01 

5. 818E+01 

10 

5.  675E+01 

5. 675E+01 

11 

5.  317E+01 

5.  317E+01 

12 

5.  585E+01 

5.  585E+01 

13 

5.  156E+01 

5. 156E+01 

14 

5.  048E+01 

5. 04BE+01 

15 

4. 744E+01 

4. 744E+01 

1C 

4. 511E+01 

4. 511E+01 

17 

4. 458E+01 

4.  458E+01 

18 

4. 314E+G1 

4.  314E+01 

19 

3. C34E4C1 

3.  634E+01 

20 

2. 6G7E+01 

2. 667E+01 

2  1 

1.  933E+01 

1.  933E+01 

22 

1.  455E+01 

1.  4  55E+01 

23 

1.  U3E+01 

1.  113E+01 

24 

8.  826E+00 

8-  82  6E+00 

25 

7.  429E+00 

7.  429E+00 

30 

2. 238E+00 

2.  238E+00 

35 

5.  890E-01 

5. 890E-01 

40 

1.  550E-01 

1.  550E-01 

45 

4.  082 E- 02 

4.  082E-02 

50 

1.  078E-02 

1.  078E-02 

70 

5.  550E-05 

5.  550E-05 

100 

1.  969E-08 

1.  969E-08 

Ar  =  1  nm  and  C[  =  0.  883  x  10"^  ,  (Figure  1).  that  is.  the  integral  over  the  whole 
size  range  becomes  1.  One  need  only  multiply  the  normalized  distribution  function 
by  the  values  of  N  in  Table  3  in  order  to  obtain  the  actual  size  distribution  at  any 
height  in  the  atmosphere. 

The  refractwe  index  for  the  aerosol  particles  has  been  assumed  real  (that  is, 
no  absorption)  for  A  s  0.  6  um,  For  A  >  0.  6  nm  the  imaginary  part  of  the  refrac¬ 
tive  index  (that  is,  absorption)  is  assumed  to  increase  linearly  to  a  value  of  0.  1 
for  A  >  2  ym.  This  modei  for  the  refractive  index  is  based  on  measurements  by 
Volz  (1957). 

The  total  numbers  (Table  3)  of  aerosol  particles  N  per  unit  volume  for  the 
"clear"  atmosphere  have  been  adjusted  so  that  the  total  extinction  coefficient  at 


i 

l 

k 


wavelength  X  =  o.  55  pm  becomes  identical  to  the  values  used  in  the  atmospheric 
attenuation  model  by  Elterman  (1988,  1970)  at  each  altitude.  This  model  corres¬ 
ponds  to  an  atmospheric  condition  of  approximately  23  km  ground  visibility.  The 
"hazy"  model  is  identical  to  the  "clear"  model  for  altitudes  II  >  5  km.  For 
H  <  5  km,  the  total  number  of  aerosol  particles  increases  exponentially  to  reach 
a  value  at  ground  level  corresponding  to  a  ground  visibility  of  approximately  5  km. 

3.  I TNIHMENTALS  OF  ATMOSPHERIC  TRANSMITTANCE 

Transmittance  of  radiation  in  the  atmosphere  is  complex  owing  to  the  de- 
penucnce  of  scattering  and  absorption  coefficients  on  a  number  of  different  physical 
properties  of  the  atmosphere  (Goody,  1964).  In  general,  we  have  for  the  mono¬ 
chromatic  transmittance  of  radiation  along  a  path  in  the  atmosphere 

t  -  exp  (-fhL|  (3) 

where  y  is  an  attenuation  coefficient  and  6  L  is  the  length  of  the  path  traversed  by 
the  radiation.  The  attenuation  coefficient  (7)  is  given  by 

y  -  a  +  k  (4) 

where  o  is  the  scattering  coefficient  and  k  is  the  absorption  coefficient.  Both  of 
these  coefficients  are  defined  in  km~*.  Equation  3  is  only  strictly  valid  when 
applied  to  monochromatic  radiation. 

The  scattering  and  absorption  coefficients  defined  in  Eq.  (4)  depend  on  both 
molecules  and  aerosols  in  the  atmospheric  path.  If  electromagnetic  radiation  is 
incident  on  a  molecule  or  an  aerosol,  a  portion  of  the  radiation  is  absorbed  and 
the  rest  is  scattered  in  all  directions.  Thus,  we  must  make  the  following  defini- 
ti.ons ' 


cr 


a 

m 


+  a 


(4a) 


k 


=  k  +  k 


(4b) 


where  the  subscripts,  m  and  a,  indicate  molecule  and  aerosol  respectively.  We 
will  endeavor  in  this  report  to  provide  quantitative  information  on  all  four  of  the 
quantities  defined  in  Eqs.  (4a)  and  (4b). 

The  molecular  scattering  coefficient  depends  only  on  the  number  density  of 
molecules  in  the  radiation  path,  whereas  the  molecular  absorption  coefficient  is 
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a  function  of  not  only  the  amount  of  absorbing  gas,  but  also  the  local  temperature 

and  pressure  of  the  gas.  The  wavelength  dependence  of  molecular  (Rayleigh) 

“4 

scattering  is  very  nearly  ~  X  .  The  variation  of  the  molecular  absorption 
coefficient  with  wavelength  is  much  more  complicated,  being  a  highly  oscillatory 
function  of  wavelength  due  to  the  presence  of  numerous  molecular  absorption  band 
complexes.  These  band  complexes  result  for  the  most  part  from  minor  atmos¬ 
pheric  constituents.  The  responsible  molecules  are  (in  order  of  importance):  H2O, 
CO.^,  O3,  N2O,  CO,  O2.  CH^  ,  N2-  All  of  these  molecules  except  H2O  and  O3 
are  assumed  to  be  uniformly  mixed  by  volume. 

The  quantities  <r  and  k  depend  on  the  number-  density  and  size  distribution  of 
aerosols  as  well  as  on  their  complex  index  of  refraction.  The  amount  of  absorbed 
radiation  and  the  distribution  of  scattered  light  can  be  derived  theoretically  for 
spherical  particles  (and  some  other  simple  forms)  by  a  theory  developed  by  Mie 
(1908).  According  to  this  theory  the  distribution  of  radiation  with  scattering  angle 
can  be  expressed  as: 


i  (</>)  =  I 
r  o 


(5) 


where  ij  and  i2  represent  the  two  polarized  components  of  the  scattered  radiation 
(van  deHulst,  1957).  They  are  functions  of  the  ratio  of  particle  size  to  wavelength, 
the  particle  refractive  index,  and  the  scattering  angle <£.  Io  is  the  intensity  of 
incident  radiation.  If  there  are  N  particles  per-  unit  scattering  volume  with  a  size 
distribution  n(r),  replacing  ij  and  i2  by  /  i-  2<r)n(r)dr  will  yield  the  integrated 
angular  scattering  intensity  \(<p)  from  the  unit  volume  of  the  particle  cloud  into  one 
steradian. 

r  4  rr 

The  integral  I  I(<£)dfl  gives  the  total  amount  cf  scattered  radiation  which 
o 

is  also  equal  to  I  '  cr  . 

1  o  a 

ICq.  (5)  can  be  written  as: 

I (4>)  =  IQ  •  a  ■  <?<<#>)  (G) 

where  v(^>)  is  the  normalized  scattering  function,  and  the  integral 

j  becomes  1. 

o 

Figure  2G  shows  ^i'p)  for  molecular  scattering  (wavelength  independent)  and 
aerosol  scattering  for  two  wavelengths  for  the  aerosol  model  described  in 
Section  2,3.  An  absorption  coefficient,  k^,  for  aerosol  particles  can  be 
defined  in  the  same  manner  as  the  scattering  coefficient,  <Ja-  The  absorption  and 
scattering  coefficients  due  to  aerosol  distributions  are  usually  smooth  functions 
of  wavelength. 
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The  molecular  absorption  coefficient  for  a  single  spectral  line  can  be  given  to 
good  accuracy  in  the  lower  50  km  of  the  atmosphere  by 


So 


it  [(v  -  rQ)"  +  o2j 


(7) 


where  S  is  the  line  intensity,  a  is  the  line  half-width,  v  is  the  frequency  of  the 
line  center  and  v  is  the  frequency  at  which  the  absorption  coefficient  is  required. 
See  the  Appendix  for  the  relationship  between  v  and  X.  The  half-width  depends  on 
the  pressure  and  temperature  approximately  as  follows: 


(8) 


The  line  intensity  also  depends  on  the  temperature  'hrough  the  Boltzmann  factor 
and  the  partition  function  (Goody,  1964  and  Penner,  1959).  Above  50  km,  Doppler 
broadening  (Goody,  1964)  should  be  considered. 

Let  us  consider  a  plane  parallel  atmosphere  and  a  frequency  at  which  several 
lines  belonging  to  one  or  more  molecular  constituents  contribute  to  the  absorption 
coefficient.  The  molecular  absorption  coefficient  in  any  one  of  the  (j)  layers  will 
be  given  by: 


k 


mj 


mjt 


(9) 


where  the  (i)  summation  is  over  a'l  molecular  absorption  lines  due  to  all  absorbing 
species  which  are  close  enough  to  the  frequency  v  to  contribute  significantly  to  the 
total  molecular'  absorption  coefficient  defined  in  Eq.  (9).  Introducing  the  molecular 
scattering  (onJ),  and  the  absorption  (kg)  and  scattering  (oQ)  due  to  aerosols  as 
defined  in  Eqs.  (4a)  and  (4b),  we  have  for  the  total  attenuation  coefficient  in  the 
(jth)  layer: 


7.  =  k  ,+ct  .  +  k  •  +  o  .  . 
.1  mj  mj  aj  aj 


(10) 


Thus,  according  to  Eq.  (3),  the  transmittance  of  an  atmosphere  having  (j)  layers 
is  given  by: 


t(v)  =  exp 


(•  4lj)' 


(U) 
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It  should  be  noted  from  Eq.  (11)  that  A L  does  not  necessarily  denote  vertical  dis¬ 
tance.  The  problem  of  slant  path  transmittance  is  easily  solved  by  writing 
ALj  =  sec  0  where  Az..  represents  the  increment  of  vertical  distance.  In 

this  case,  Eq.  (11)  becomes 

tm0/)  ■  [r(r)]M  ('2) 

where  M  =  sec  0  is  the  secant  of  the  angle  between  the  vertical  and  the  direction  of 
observation.  Thus  it  becomes  very  easy  to  convert  vertical  transmittance  values 
to  slant  path  transmittance.  As  9  approaches  90°,  the  result  becomes  meaningless 
due  to  the  plane  parallel  atmosphere  assumption  (see  Section  7). 

All  of  the  above  discussion  is  valid  only  at  a  single  wavelength  (or  very  nar¬ 
row  wavelength  interval  which  is  known  to  be  small  compared  with  the  width  of  any 
absorption  lines  in  the' vicinity).  A  sufficiently  narrow  spectral  interval  for  most 
purposes  is  0.001  wavenumber  (cm'*).  The  process  which  so  limits  the  validity 
of  the  above  relations  is  molecular  absorption.  As  noted  above,  krn(i')  is  a  highly 
oscillatory  function  of  frequency  whereas  k  (v),  a  (v),  and  a  (v)  are  much  more 
slowly  varying  functions- 

When  it  is  desired  to  calculate  the  atmospheric  transmittance  over  a  large 
spectral  interval,  it  is  necessary  to  average  Eq.  (11)  or  Eq.  (12)  over  the  desired 
spectral  interval  (see  Eq.  (I4)j .  When  many  absorption  lines  are  present,  this 
can  be  a  particularly  tedious  job  (even  on  a  computer).  Various  techniques  have 
been  developed  to  treat  this  problem  and  are  discussed  in  Section  5. 

4.  ATMOSPHERIC  TRANSMITTANCE  AT  SELECTED  FREQUENCIES 

4.1  List  of  Laser  Frequencies 

It  is  necessary  to  limit  the  number  of  specific  frequencies  at  which  we  quote 
atmospheric  absorption  and  scattering  coefficients.  We  have  therefore  chosen  the 
list  of  specific  laser  wavelengths  (Table  4)  ai  which  values  of  these  coefficients 
are  given  in  Table  5. 

4.2  Scattering  and  Absorption  Coefficients 

Each  atmospheric  model  defined  in  Section  2  has  been  used  as  a  basis  for 
calculating  the  absorption  and  scattering  coefficients  given  in  Table  5  for  each 
laser  wavelength.  Numerical  values  for  the  coefficients  were  computed  for  the  0, 

1,  2,  etc.  ,  km  levels  and  then  logarithmically  interpolated  for  each  atmospheric 
layer.  The  aerosol  attenuation  coefficients  given  in  Table  5,  which  correspond  to 
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Table  4.  Laser  Wavelengths-' 


0.  33V  1  pm 

Nitrogen 

0.  4880  pm 

Argon 

0.  5145  pm 

Argon 

0.  6328  pm 

Helium  -Neon 

0.  6943  pm 

Ruby 

0.  86  pm 

Gallium  Arsenide 

1.06  pm 

Neodymium  in  GLass 

1.  536  pm 

Erbium  in  Glass 

3.  39225  pm 

Helium  Neon 

10.591  pm 

Carbon  Dioxide- 

27.  9  pm 

Water  Vapor 

337  pm 

Hydrogen  Cyanide 

the  clear  and  hazy  atmospheres  defined  in  Section  2.3,  are  identical  for  all  geo¬ 
graphical  and  seasonal  models.  In  Table  5,  the  four  extinction  coefficients  are 
listed  as  a  function  of  altitude  for  each  model:  km  and  om  are  the  molecular  ab¬ 
sorption  and  scattering  coefficients,  and  kg  and  oa  are  the  aerosol  absorption  and 
scattering  coefficients.  If  the  total  attenuation  coefficient  is  required,  the  four 
coefficients  can  be  summed: 


+  a 


+  k  +  a 
a  a 


(13) 


4.3  Interpolation  Procedures 

In  general,  it  is  permissible  to  interpolate  linearly  with  respect  to  wavelength 
between  tabulated  values  in  Table  5  except  for  molecular  absorption  coefficient 
(k^).  This  is  true  because  of  the  smooth  variation  of  scattering  coefficients  and 
of  aerosol  absorption  coefficients  with  wavelength,  and  the  rapid  spectral  variations 
in  molecular  absorption  coefficients.  The  tables  have  been  constructed  so  that 
linear  interpolation  with  respect  to  height  is  expected  to  give  correct  values  to 
±8  percent.  If  a  given  atmospheric  condition  differs  from  one  of  the  ten  models  of 
Section  2,  larger  errors  may  result. 

Again  it  should  be  remembered  that  one  should  never  attempt  to  interpolate 
with  respect  to  wavelength  between  molecular  absorption  coefficient  values  (k^) 
given  in  Table  5, 

*  These  lasers  may  emit  over  a  range  of  wavelengths.  However,  the  computed 
values  in  Table  5  were  derived  for  the  specific  wavelength  values  of  Tabic  4. 
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Table  5  (Contd).  Values  o'?  Attenuation  Coeffic tent/km  £s  a  Function  of  Altitude  for  Each  laser  Wavelength  and 
Each  Atmospheric  Model  in  Section  2  (k^  =  molecular  a.bsorption.  =  molecular  scattering  kg  =  aerosol 
absorption,  <r  -  aerosol  scattering) 
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Each  Atmospheric  Model  in  Section  2  (km  =  molecular  absorption,  =  molecular  scattering,  =  aerosol 
absorption,  a,  =  aerosol  scattering) 
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Having  obtained  the  attenuation  coefficients  y[Eq.  (13)],  the  laser  transmit¬ 
tance  for  a  horizontal  path  of  length,  L,  at  height,  h,  is  given  by  t  =  exp  (-7.L) 
when  the  suffix  j,  refers  to  the  layer  at  height,  h.  For  slant  paths  between  two 
altitudes,  Eq.  (11)  should  be  used  with  the  appropriate  correction  for  zenith  angle. 


5,  ATMOSPHERIC  TRANSMITTANCE  AT  LOR  RESOLUTION  FROM 
0.25  TO  23 


The  calculations  presented  in  Section  4  were  made  at  infinite  resolution 
(monochromatic  calculations).  In  practice  it  is  impossible  to  measure  radiation 
transmitted  at  a  single  frequency.  Instead  one  measures  the  transmittance  r.  v(v) 
averaged  over  the  spectral  interval,  Au  ,  accepted  by  the  receiver,  as  indicated 
in  the  equation 


wv) 


(v)  dr, 


(14) 


where  r  is  the  central  frequency  in  the  interval,  Av.  Consequently  for  many  appli¬ 
cations  one  is  interested  in  knowing  the  transmittance  of  the  atmosphere  averaged 
over  a  relatively  wide  spectral  interval,  that  is,  for  low  resolution. 

Thus  the  term  transmittance  is  somewhat  ambiguous  unless  it  is  qualified  by 
some  indication  of  the  spectral  resolution,  At1,  over  which  it  is  averaged.'  This  is 
particularly  true  in  the  case  of  molecular  absorption,  since  the  absorption  coeffi¬ 
cient  kJn  is  a  rapidly  varying  function  of  frequency.  It  is  because  of  the  rapid 
variation  of  kjj,  with  frequency  that  the  averaged  transmittance  r  does  not,  in 
general,  obey  the  simple  exponential  law  given  in  Eq.  (3).  That  is. 


T&v<v) 


dr  *  exp  (-  K(v)AL) 


(15) 


where  represents  the  net  monochromatic  molecular  absorption  coefficient  as 
given  in  Eq.  (9)  and  where  K  is  an  average  absorption  coefficient,  which  cannot  in 
most  cases  be  defined  when  Ar  is  much  greater  tlian  the  half  width  of  a  spectral 
line  [see  Eq.  (8)] . 

On  the  other  hand  the  molecular  scattering  coefficient  (ari)  and  the  aerosol 
scattering  and  absorption  coefficients  [eta  and  1^,  see  Eqs.  (4a)  and  (4b)j  are  slow¬ 
ly  varying  functions  of  frequency,  and  the  average  transmittance  (which  can  easily 
be  obtained  by  interpolation  from  the  monochromatic  values,  see  Figures  21  and  22) 
obeys  the  simple  exponential  law  |Eq.  (3)J  provided  only  the  direct  transmitted 
beam  is  being  observed  (see  Zuev  et  al,  1967). 


I 
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There  are  four  basic  approaches  to  obtaining  low  resolution  transmittance 
values  for  a  given  path  through  the  atmosphere  due  to  molecular  absorption.  These 
are:  (1)  direct  measurements  over  the  required  path,  (2)  measurements  in  the 
laboratory  under  simulated  conditions,  (3)  line-by-line  (monochromatic)  calcula¬ 
tions  based  on  detailed  knowledge  of  spectroscopic  line  parameters  (see  Section  3) 
which  are  then  averaged  over  the  required  spectral  interval,  and  (4)  calculations 
based  on  band  modeL  techniques  (which  use  available  laboratory  and/or  field  trans¬ 
mittance  measurements  or  actual  line  data  as  a  basis). 

From  the  point  of  view  of  computations,  method  3  involves  a  considerable 
amount  of  work  and  computer  time,  and  consequently  method  4  has  been  used  most 
frequently. 

In  this  section  we  will  briefly  discuss  the  problem  of  continuum  absorption  in 
the  atmospheric  window  region  (Section  5.  1).  Then  we  will  summarize  the  appli¬ 
cation  of  band  models  to  molecular  absorption  calculations  (section  5.2)  and 
present  a  detailed  description  of  a  graphical  prediction  scheme  for  determining 
the  transmittance  of  the  atmosphere  due  to  line  absorption,  continuum  absorption, 
and  scattering  (average  over  20  cm"*)  for  a  given  atmospheric  path.  First,  the 
reader  is  provided  with  figures  for  determining  the  relevant  absorber  amounts  in 
the  required  atmospheric  path  (Section  5.  3),  which  are  then  used  in  conjunction 
with  the  transmittance  charts  (Section  5.4)  to  determine  (separately)  the  transmit¬ 
tance  due  to:  (1)  molecular  (line)  absorption  (H^O,  O3,  end  the  umlorrniy  mixed 
bases).  (2)  molecular  (continuum)  absorption  (H20,  N2),  (3)  molecular-  (Rayleigh) 
scattering,  and  (4)  aerosol  extinction  (absorption  and  scattering),  as  a  function  of 
frequency  (wavelength). 

Finally  the  reader  has  to  multiply  together  the  individual  transmittances  at 
each  wavelength  in  order  to  determine  the  total  transmittance  at  that  wavelength. 
That  is; 

v^y(v)  (Total)  =  T^y(i/)  (Une  absorption)  x  t  (1/)  (continuum  absorption) 

X  (Rayleigh)  x  *£*,<«')  (aerosol)  . 

5.1  Continuum  Absorption  in  the  Atmospheric  Window  Regions 

There  are  a  number  of  spectral  regions  of  weak  or  negligible  absorption  by 
atmospheric  molecules  in  the  infrared  which  are  generally  called  atmospheric 
windows.  Various  attempts  have  been  made  to  determine  the  transparency  of  these 
window  regions  (Yates  and  Taylor,  1960;  and  Streete,  1968).  Yates  and  Taylor 
(1960)  concluded  that  the  opacity  in  the  near  infrared  windows  is  limited  at  sea 
level  by  aerosol  scattering.  This  conclusion  is  consistent  with  the  results  which 
are  derived  from  our  aerosol  chart  presented  in  Section  5.4 
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As  there  are  no  laboratory  data  presently  available  from  which  the  molecular 
absorption  between  absorption  bands  can  be  inferred  (with  two  important  excep¬ 
tions),  we  will  here  assume  that  the  gaps  in  the  charts  presented  in  Figures  IB 
through  18  are  real  gaps  ol  zero  molecular  absorption.  Thus,  the  attenuation  in 
these  regions  is  dependent  entirely  on  molecular  and  aerosol  scattering. 

However,  laboratory  absorption  data  are  available  for  the  water  vapor  con¬ 
tinuum  (Burch,  1970;  and  McCoy  et  al,  1969).  from  8-13pm  and  also  for  the  pres¬ 
sure  induced  (continuous  absorption)  band  of  nitrogen  (Reddy  and  Cho.  1965) 
centered  near  4.  3  The  strongest  part  of  this  nitrogen  band  is  not  too  import¬ 

ant  as  it  is  centered  beneath  the  very  strong  4.  3  pm  CQ2  band.  However,  the 
nitrogen  absorption  extends  out  beyond  the  region  of  C02  absorption  near  4.  0  um 
and  becomes  predominant,  accounting  for  2  0  percent  absorption  in  a  10  km  path 
at  sea  level. 

Due  to  the  continuous  nature  of  these  absorbing  regions,  the  transmittance 
follows  a  simple  exponential  law  which  is  represented  by  the  scale  in  Figures  19 
and  20. 


5,2  Outline  of  Uund  Model  Techniques 

Basically,  a  band  model  assumes  an  a -ray  of  lines  having  chosen  intensities, 
half-widths  and  snacings,  winch  can  be  adjusted  to  represent  the  line  si  rue  tut  e  in 
some  part  of  a  real  band.  For  a  particular  banc!  model  the  mean  transmittance  can 
be  represented  by  a  mathematical  expression  (transmittance  function);  for  example 


expressed  in  terms  of  pressure,  P,  effective  path  length  (or  absorber  concentra¬ 
tion)  AL,  and  one  or  more  frequency  dependent  absorption  coefficients  (C(i')  s). 

Several  band  models  have  been  developed  (see  Goody,  1952  and  1964;  Elsasser, 
1938;  Plass,  1958;  King,  1959  and  1964;  and  Wyatt  et  al,  1964);  the  Elsasser  model 
(1942)  and  the  Goody  model  (1952  and  1964)  being  moot  well  known.  The  Elsasser 
band  model  consists  of  an  array  of  equally  spaced  Identical  lines  of  Lorentz  shape 
(see  Eq.  (7)j  .  This  model  has  been  applied  to  absorption  bands  which  have  a 
regular  line  structure,  for  example,  to  some  bands  of  C02,  N20,  CO,  CHq  and 
O,  .  The  Goody  model,  or  the  other  hand,  assumes  that  the  band  is  composed  of 
spectral  lines  with  an  exponential  intensity  distribution  and  with  random  spacing 
between  lines.  Again  the  lines  are  assumed  to  have  a  Lorentz  shape.  This  model 
has  generally  been  applied  to  bands  which  have  an  ii  regulai  line  structure,  for 

example,  to  some  H,pO  and  bands. 

In  practice  the  wavelength  dependent  coefficients  are  determined  for  each 
absorbing  gas  separately  using  laboratory  transmittance  values  measured  under 
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known  conditions,  that  is  by  solving  Eq.  (16)  for  the  CiD's  at  each  v  for  known 
values  of  t(i/)(  AL  and  P. 

The  absorption  coefficients  obtained  in  this  way  are  then  used  in  the  band 
model  transmittance  function  [Eq.  (16))  to  determine  t(i/)  as  a  function  of  frequency 
v  for  other  values  of  AL  and  P  for  each  absorber.  Finally,  the  total  mean  trans¬ 
mittance  for  molecular  absorption  is  given  by  the  product  of  the  mean  transrnit- 
tances  of  the  individual  absorbers  at  each  frequency. 

t  (v)  (molecular  absorption)  *  t  (v  )  (11^0)  x  t^(i/)  (CO.,) 


x  T^(N2°'^Ar(,,>  <°3> 


(17) 


Exact  analytical  expressions  have  been  obtained  for  most  of  the  band  models. 
However,  they  are  sometimes  difficult  to  use  and  simpler  approximations  have 
been  found  to  apply  in  two  limiting  conditions,  which  are  common  to  all  band 
models.  These  simpler  expressions  are  the  well  known  "weak  line"  and  "strong 
line"  approximations  (see  Goody,  1964;  and  Plass,  1968)  for  which  the  transmit¬ 
tance  is  a  function  of  the  absorber  amount.  Ah,  and  the  product  of  the  pressure  and 
absorber  amount,  PAL,  respectively  (for  a  given  temperature). 

The  weak  line  approximation,  which  corresponds  to  the  exponential  law,  is 
valid  when  the  absorption  is  small  at  the  line  centers  (generally  for  high  pressures 
and  low  absorber  amounts).  Unfortunately  this  case  is  rarely  applicable  to  condi¬ 
tions  existing  in  the  terrestrial  atmosphere. 

The  strong  line  approximation  is  applicable  where  the  lines  are  completely 
absorbing  at  their  centers;  the  effect  of  increasing  the  amount  of  absorber  is  then 
confined  to  the  edges  or  wings  of  the  lines. 

The  regions  of  validity  of  the  strong  and  weak  line  approximations  for  the 
Elsasser  and  Goody  models  are  discussed  in  Plass  (1958). 

For  practical  purposes  most  problems  fall  in  either  the  strong  line  approxima¬ 
tion  region  or  the  intermediate  regime.  It  is  perhaps  interesting  to  compare  the 
functional  torm  of  the  strong  line  approximation  for  some  of  the  modeLs. 

For  the  Elsasser  Model 


7&v{v)  "  1  '  erf  AL.PJ 


where  erf  is  the  error  function 


1/2 


erf  (z) 


.  JL  f 


e'1*  dt 


(18) 


and  C(r)  is  an  average 


•*  z 

rtiorption  coefficient  which  can  be  written  ir.  terms  of  average  values  of  the  mole¬ 
cular  line  constants  detined  in  Section  3. 
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The  strong  line  approximation  for  the  Goody  model  is  given  by 

=  exp  {(C'(v)  ALP]1/2}  (19) 

where  C'(r)  is  defined  in  a  similar  manner  to  C(r)  above. 

King  (1959)  attempted  to  generalize  the  above  expressions  by  writing 

(20) 

so  that  when  n  is  set  equal  to  zero  cr  unity,  the  average  transmittance  agrees 
with  the  weak  or  strong  line  approximations  respectively. 

In  this  report  we  have  determined  from  laboratory  and  synthetic  transmittance 
data,’;-  three  empirical  transmittance  functions  based  on  Kq.  (20)  for  (1)  H2O. 

(2)  Og.  and  (3)  the  combined  contributions  of  the  uniformly  mixed  gases.  It  was 
found  that  Kq.  (20)  gave  the  best  fit  to  laboratory  and  theoretical  data  over  a  wide 
spectral  interval  when:  (1)  n  =  0.  9  for  H.,0,  (2)  11  =  0.  4  for  O3,  and  (3)  n  *  0.75 
for  the  uniformly  mixed  gases. 

The  corresponding  empirical  transmittance  functions  were  found  to  give  better 
agreement  with  laboratory  and  synthetic  transmittance  data,  than  the  commonly 
used  band  models  over  a  wide  range  of  pressures  and  absorber  amounts. 

The  procedure  used  for  determining  the  parameter  »  and  function,  f,  defined 
in  Kq.  (20)  will  be  briefly  outlined  below.  By  taking  the  logarithm  of  the  inverse 
of  Kq.  (20),  it  will  be  seen  that 

n  log  P  +  log  aL  =  logy  1[T£I/(r)j  -  log  C"  (r).  (21) 

Thus  for  a  given  frequency  v  and  fixed  values  of  average  transmittance 
t^i/(u),  the  righthand  side  of  Kq.  (20)  becomes  a  constant. 

A  mean  value  for  n  was  determined  from  Kq.  (21)  for  a  wide  range  of  fre¬ 
quencies  and  several  values  of  T.  Then  for  each  frequency,  T  was  plotted  against 
log  AL  Pn,  and  the  curves  superimposed  and  the  best  mean  curve  determined. 

The  "mean"  curve  thus  obtained  constitutes  the  empirical  transmittance  function 
and  is  displayed  as  the  transmittance  scale  and  associated  scaling  factor, 
log  (AKPn),  in  the  prediction  cha-ts  given  in  Figures  10  through  18. 

So  far  the  discussion  has  been  restricted  to  constant  pressure  paths.  In  the 
following  section  we  will  discuss  the  application  of  band  models  to  variable  pressure 
paths. 

♦Synthetic  spectra  her  e  refers  to  monochromatic  (line  by  line)  transmittance 
calculations  degraded  in  resolution  to  20  cm”1. 
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5.3  I', qui  valent  C.onslant  PresMiro  Paths  und  Absorber  Amounts 

Atmospheric  paths  may  be  divided  into  two  types,  namely,  constant  pressure 
paths  (that  is,  horizontal  paths  of  up  to  about  100  km)  and  slant  or  variable  pres¬ 
sure  paths.  The  extent  of  the  former  types  of  path  is  limited  due  to  the  curvature 
of  the  earth.  Since  the  transmittance  functions  for  constant  pressure  paths  arc 
Simpler  to  handle  than  those  for  variable  pressure  paths,  the  problem  is  simplified 
by  reducing  the  slant  path  to  an  "equivalent  constant  pressure  path";  that  is,  to  a 
constant  pressure  path  which,  ideally,  has  the  same  transmittance  at  each  wave¬ 
length  as  does  the  actual  slant  path. 

In  tliis  section  wc  will  discuss  the  construction  and  application  of  charts  for 
determining  the  equivalent  sea  level  path  quantities  for  molecular  absorption  (in¬ 
cluding  continuum  absorption  in  the  4  pm  and  8-13  pm  region),  molecular  scatter¬ 
ing  and  aerosol  extinction  based  on  the  model  atmospheres  presented  in  Section  2. 

5.3.1  MOLKCUI.AR  ABSORPTION 


It  has  been  assumed  that  for  most  practical  purposes,  the  average  trans¬ 
mittance  is  dependent  on  the  product  Al.Pn  where  AL  is  the  absorber  amount,  and 
P  is  the  pressure. 

Here  we  will  define  the  equivalent  sea  level  amounts  for  each  of  the  five  model 
atmospheres  given  in  Section  2  as  the  value  of  the  above  product  (reduced  to  STP) 
averaged  over  a  given  atmospheric  path  for  a  particular  absorber,  that  is 
/  P  \n 

A  Pp — j  where  the  subscript  (o)  refers  to  STP. 


Since  liyO  and  Oy  are  distributed  nonuniformly  in  the  atmosphere,  we  must 
consider  these  gases  separately  from  the  uniformly  mixed  gases. 

For  water  vapor,  n  is  set  equal  to  0.0  and  AL0  *  Kw  where 

w(z)  is  the  water  vapor  concentration  at  altitude,  z,  in  units  of  gm  cm"2/lun 
(obtained  from  Table  1),  and  rt  is  the  range  in  km.  Since  w(z)  is  giver  in  mass 
units  there  is  no  temperature  or  pressure  correction  necessary  to  reduce  this 


quantity  to  STP.  The  variation  with  altitude  of  w 


/pX0.9  /•*  /p\0.9 

\Fo/  "'HO 


cb 


(denoted  by  wh  and  wv  respectively)  is  shown  in  Figures  2  and  3.  The  subscripts, 
h  and  v,  refer  to  horizontal  and  vertical  paths  respectively. 

For  the  uniformly  mixed  gases,  n  is  set  equal  to  0.  75,  but  the  absorber  con¬ 
centration,  AL0  (reduced  to  STP)  for  a  horizontal  path  of  range  B  at  altitude  z 

is  given  by  AL.0  =  cR  m  where  c  is  the  fractional  concentration  by  volume  of 

the  absorbing  gas  (see  Table  2  and  Appendix).  The  terms  P  and  T  refer  to  the 
pressure  and  temperature  at  altitude  z.  Therefore,  we  have  in  this  case 
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AL0(p~)  (t)  ’  The  variation  with  altitude  of  (~^j  and 

f°°  /  p\3-75  /T0\ 

J  YT/  Z  (dcnoted  by  <h  and  |y  respectively), is  shown  in  Figures 


4  and  5. 


.0.4 


For  ozone,  n  is  set  equal  to  0.4  (see  Section  5. 2),  andAL0^^-J*  =  u(z) 

where  u(z)  is  the  ozone  concentration  in  units  of  (cm  -  atm)/km  at  STP  (obtained 

/pip.4  /•«  /P\0-4 

from  Table  1).  The  variation  with  altitude  of  u(z) ^5- j  and  J  u(z)^p-^  dz 

(denoted  by  u^  and  uy  respectively)  is  shown  in  Figures  6  and  7. 

5.  3.  2  NI  TROGEN  CONTINUUM  (4  pm  REGION) 

The  absorption  due  to  the  nitrogen  collision  induced  band  is  proportional  to 
the  square  of  the  pressure.  The  effective  absorber  amount  (reduced  to  STP)  for 


a  horizontal  path  of  range  R  at  altitude  z  is  cR 


m 


The  variation  with 


,  p  x2  (T  \  IV  \2  /T  \ 

*  "  1  ’  0  ■■'■■■  !  '  -  ■‘ttl  paths  and  /  (p~)  {-=?!  dz  for  vertical  paths 

•At  Vo/  x  1  ' 


I  IOI  llUi  UOludl 


altitude  ot^J  yT/ 

(denoted  by  n^  and  nv  respectively)  is  shown  in  Figures  8  and  9. 


5.  3.  3  WATER  VAPOR  CONTINUUM 

Laboratory  measurements  of  the  water  vapor  continuum  indicate  that  self¬ 
broadening  due  to  water  vapor  molecules  is  much  niore  important  than  broadening 
due  to  other  atmosphei  ,c  molecules  (Burch,  1970).  The  appropriate  equivalent 
sea  level  absorber  amount  for  the  water  vapor  continuum  is 


bh  =  P 


Pp  O 
L  K2U 


+  0.  005  It 


\r 


Y1  .  > 

Ph2o)J  w'*; 


for  horizontal  paths,  where  P  is  the  total  pressure  and  Pj^O  *s  the  Par,-'a^  pressure 
of  water  vapor  (Pp^O  =  4-'’bxl0  S  w(z)T(z)  atm,  where  th.e  constant  is  given  by  the 
ratio  of  the  universal  gas  constant  to  the  molecular  weight  of  water  vapor)  and  w(z) 
is  the  number  of  (gm  crn'^/km)  of  water  vapor  at  altitude  z  obtained  from  Table  1. 
The  equivalent  sea  level  amount  for  vertical  paths  is  given  by 


r®  1 

bv(z)  J  ~P 


z  Po 


11,0 


+  0.005 


(P‘P'V>)' 


w(z)  dz 


The  variation  with  altitude  of  these  functions,  denoted  by  b^  and  by  respectively, 
is  shown  in  Figures  10  and  11. 

5.3.4  MOLECULAR  SCATTERING 


Molecular  scattering  is  a  function  of  atmospheric  density  p( z)  (Section  3). 

For  each  of  the  five  model  atmospheres  given  in  Section  2,  Figures  12  and  13  show 

/•“ 

the  variation  with  altitude  z  of  p(z)/p  and  /  p(z )/ pQ  dz  (denoted  by  and  mv 

~4  v 

respectively)  where  pQ  is  the  density  of  air  at  STP. 

5.3.5  AEROSOL  EXTINCTION 


Aerosol  extinction  is  a  function  of  the  height  variation  of  the  aerosol  number 
density  N(z)  for  a  given  size  distribution  function  (Section  2.  3).  Figures  14  and  15 


show  the  variation  with  altitude  of 


N(z) 


and 


/N(z)dz 

V 


(denoted  by  a,  and  a 
h  v 


respectively)  for  the  two  hazy  models  given  in  Section  ?,.  3  where  N0  is  the  number 
density  (cm^)  of  aerosols  at  sea  level  for  a  visibility  of  23  km. 


5.  3.  6  UTILIZATION  OF  EQUIVA-LENT  PATH  FIGURES 


For  a  given  atmosphe*  ic  path,  let  us  consider  how  to  obtain  the  equivalent 
sea  level  quantities  for  HgO.Og,  the  uniformly  mixed  gases,  the  nitrogen  and 
waiti  coiitinua,  molecular  scattering  and  aerosol  extinction  which  will  be  denoted 
by  W,  U,  L,  N,  B,  M,  and  A  respectively.  At  this  point  the  reader  need  not  be 
concerned  with  the  units  of  these  quantities  which  were  chosen  to  facilitate  the  use 
of  the  prediction  charts  given  in  Section  5.4.  The  subject  of  units  will  be  discussed 
separately  in  the  Appendix. 

Figures  2  through  15  show  the  variation  with  altitude  of  the  relevant  quantities 
for  both  horizontal  and  vertical  paths  (distinguished  by  suffixes  h  and  v  respec¬ 
tively)  The  use  of  these  figures  can  best  be  illustrated  by  the  following  examples: 

5.3.6.  1  Horizontal  Paths 


To  determine  the  equivalent  sea  level  quantities  for  a  horizontal  path  of 
length  R  (<  100  km)  at  an  altitude  z,  simply  multiply  each  of  the  quantities  w^, 
uh‘  ”h‘  bh'  mh  ant*  ah  ^b^ined  from  Figures  2,  4,  6,  8,  10,  12,  14  for  the  ap¬ 
propriate  altitude  z),  by  R  (in  kms).  For  example  in  the  case  of  molecular  line 
absorption  by  water  vapor,  W  is  given  by 


5.  3.  6.  2  Slant  Paths 


(1)  The  equivalent  sea  level  quantities  (W ,  L,  U,  N,  B,  M  a  d  A)  for  a  slant 
path  from  altitude  z  to  the  top  of  the  atmosphere  at  zenith  angle  are  determined 
by  multiplying  each  of  the  quantities  wv>  <v>  uv,  n  ,  bv<  mv>  and  a^  (obtained 
from  Figures  3,  S,  7,  9,  11,  13,  15  for  the  appropriate  altitude  z)  by  sec  0  (where 
0  <  80°).  For  example,  in  the  case  of  molecular  line  absorption  by  water  vapor 

W  =  wy(z)  sec  0  .  (22) 

For  0  >  80°,  sec  0  should  be  replaced  by  M(0)  given  in  Figure  32  (see  Section  7). 

(2)  For  a  slant  path  of  range  R  between  altitudes  z ,  and  z„,  the  value  of  W  for 
molecular  line  absorption  by  water  vapor  is  given  by 

W  =  |wv(Zl)  -  wv  <z,)j  x  -  .  (23) 

Similar  expressions  hold  for  I,,  U,  N,  B,  M  and  A. 

Corrections  can  also  be  applied  to  account  for  atmospheric  refraction  and 
earth  curvature  and  these  will  be  briefly  summarized  in  Section  7. 

Having  now  determined  W  (gm.  cm"'*),  U  (cm-atm  at  STP),  L,  (km),  N  (km), 

B  (gm  cm'2),  M  (km)  and  A  (km)  one  can  use  tne  prediction  ciiaris  given  in 
Figures  16  through  22  (Section  5.4). 

5.1  Description  and  Use  of  Prediction  (.Kurts 

The  purpose  of  this  section  is  to  present  a  simple  method  of  predicting  atmos¬ 
pheric  transmittance  (at  low  resolution)  which  is  applicable  over  a  wide  spectral 
interval  and  for  a  wide  range  of  atmospheric  paths.  Of  the  methods  available,  the 
one  originally  devised  by  Altshuler  (1961)  is  the  most  comprehensive  that  meets 
the  above  requirements,  although  based  on  fairly  old  experimental  data.  In  view 
of  this,  an  attempt  has  been  made  here  to  further  simplify  Altshuler's  method  and 
to  update  it  by  (1)  using  more  recent  laboratory  and  field  transmittance  data  (Burch 
et  al,  1962,  1964,  1965a,  1965b,  1965c,  1965d,  1967,  and  1968)  and  (2)  making 
some  line  by  line  calculations  for  molecular  absorption  (as  discussed  in  Section 
5.4)  and  averaging  over  the  spectral  interval  (A v  =  20  cm''). 

Charts  for  calculating  the  transmittance  due  to  molecular  line  absorption  by 
water  vapor  (I^O).  ozone  (O,^),  the  uniformly  mixed  gases  (COj-  NgO,  CH^,  CO, 
Og),  the  nitrogen  4  pm  continuum,  the  water  vapor  8-13  pm  continuum,  molecular 
scattering,  and  aerosol  extinction  and  absorption  are  presented  in  Figures  16 
through  22.  It  will  be  noted  that  the  effects  of  the  different  uniformly  mixed  gases 
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have  been  combined  into  one  chart  (Figure  17)  (assuming  the  mixing  ratios  given 
in  Table  2). 

A  description  and  outline  of  the  use  of  the  charts  will  now  be  given. 

Each  figure  shows  the  variation  of  an  attenuation  coefficient  with  respect  tc 
wavelength,  relative  to  a  horizontal  datum  line  for  a  given  effective  sea  level  path 
length  or  absorber  amount.  Also  Ghown  in  each  figure  is  a  vertical  transmittance 
scale  associated  with  a  set  of  scaling  factors  (for  the  absorber  amount).  In  order 
to  use  these  charts  efficiently,  it  is  necessary  to  trace  the  transmittance  scale  and 
associated  scaling  factors  for  each  of  the  figures  on  to  some  transparent  paper. 

If  the  transmittance  scale  is  now  placed  on  its  appropriate  chart  so  that  the 
scaling  factor  1  rests  on  the  datum  line  (as  shown  in  Figure  16(b)),  then  one  is  in  a 
position  to  read  off  the  transmittance  (for  1  gm.  cm"^  of  H2O)  as  a  function  of 
wavelength.  Simply  move  the  chart  horizontally  until  the  transmittance  scale  is 
vertically  above  the  required  wavelength  and  read  off  the  transmittance  where  the 
scale  crosses  the  curve.  For  example,  in  Figure  16(b)  the  scale  is  set  at  21 
microns  and  indicates  a  transmittance  of  0.  10  for  W  =  1  gm  cm'2  of  HgO.  *  If  the 
transmittance  is  required  for  a  different  amount  of  absorber,  for  example, 

W  =  10  gm  cra'^,  just  displace  the  transmittance  scale  vertically  until  the  scaling 
factor  10  is  coincident  with  the  datum  line.  Then  move  the  scale  horizontally  to 
the  required  wavelength  and  read  off  the  transmittance  as  described  above. 

The  same  procedure  applies  to  Figures  17  through  22,  25  and  33,  where  it 
will  be  noted  the  datum  line  for  scaling  the  absorber  amounts  is  expressed  as 
1  (cm-atTn)(,Tp  for  (Figures  18  and  25),  1  gm  cm'2  for  the  water  vapor  con¬ 
tinuum  (Figure  20),  and  1  km  for  each  of  the  other  charts. 

Thus  by  displacing  the  transmittance  scale  horizontally  according  to  wave¬ 
length  and  vertically  according  to  the  appropriate  equivalent  sea  level  quantity  (W, 

U,  L,  N,  B,  M  or  A)  for  Figures  16  through  22  respectively,  one  can  construct  a 
low  resolution  transmittance  spectrum  (for  Aw  2:20  cm'*). 

A  method  for  calculating  the  equivalent  sea  level  absorber  amounts  W,  U,  L, 

N,  B,  M  and  A  respectively  (in  the  appropriate  units)  for  a  given  atmospheric  path 
has  been  described  in  Section  0.  3  uaSc  d  on  the  ter.  model  atmospheres  defined  in 
Section  2.  In  the  case  where  the  relevant  meteorological  parameters  are  known 
for  the  atmospheric  path,  these  quantities  can  be  calculated  by  an  alternative 
method  described  in  the  Appendix. 

5.5  Absorption  in  the  Region  0. 25  to  0.75  pm 

The  only  atmospheric  molecule  considered  In  this  report  in  the  spectral  region 
from  0.25  to  0.70pm  is  ozone.  Weak  absorption  bands  due  to  water  vapor  and  oxygen 

♦  This  corresponds  to  the  amount  of  water  vapor  in  a  )  krn  path  at  sea  level  for 
T  =  15°C  and  78  percent  relative  humidity  (see  Appendix). 
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also  appear  in  this  region  but  because  of  a  lack  of  reliable  data  have  not  been  in¬ 
cluded  here.  Since  the  absorption  coefficient  due  to  ozone  in  this  wavelength  range 
is  independent  of  pressure,  it  is  sufficient  to  define  the  amount  of  ozone  in  an 
atmospheric  path  with  a  single  parameter  proportional  to  the  total  number  of  mole¬ 
cules  in  the  path.  Figure  23  gives  the  ozone  amount  as  a  function  of  altitude  in 
<cm-atm)g^,p/km  for  each  of  the  five  atmospheric  models,  and  should  be  used  to 
determine  the  amount  of  ozone  in  a  horizontal  pa>-h.  Figure  24  provides  total  ozone 
amounts  between  the  top  of  the  atmosphere  and  the  indicated  altitudes  for  each 
model  atmosphere.  Ozone  amounts  for  vertical  layers  within  the  atmosphere  can 
be  found  by  subtraction  of  vaiueg  read  off  Figure  24  for  two  different  altitudes. 

Slant  path  amounts  are  found  by  multiplication  of  the  vertical  amount  by  sec  ©  or 
MW)  (see  Section  7).  For  a  discussion  of  the  unit,  (cm-atm)g^,p,  used  in  this 
section,  see  the  Appendix. 

The  amount  of  ozone  in  the  required  atmospheric  path  is  determined  from 
Figures  23  and  24.  Then  the  atmospheric  transmittance  through  the  path  is  found 
as  a  function  of  wavelength  from  Figure  25  obtained  from  Inn  and  Tanaka  (1953). 

The  application  of  this  figure  is  identical  to  that  described  earlier  for  the  prediction 
charts  given  in  Section  5.4. 

S.6  Accuracy  of  Low  Resolution  Transmittance  Calculations 

The  overall  accuracy  in  transmittance  which  this  technique  provides  is  better 
than  10  percent.  The  largest  errors  may  occur  in  the  distant  wing6  of  rather 
strongly  absorbing  bands.  The  reason  for  this  error  is  that  a  unique  spectral 
curve  as  given  in  Figures  10-25  cannot  be  defined  for  all  possible  atmospheric 
paths.  The  curves  presented  in  these  figures  were  constructed  for  moderate 
atmospheric  paths  and  will  tend  to  overestimate  the  transmittance  for  very  long 
paths  and  underestimate  the  transmittance  for  very  short  paths- 

As  the  transmittance  approaches  1.  0,  the  percent  error  in  transmittance 
decreases  toward  zero,  but  the  uncertainty  in  the  absorptance  (or  emissivity) 
increases  greatly.  In  fact,  if  the  computed  transmittance  exceeds  0.96,  the  un¬ 
certainty  in  absorptance  becomes  so  great  that  the  technique  should  not  be  used 
for  a  determination  of  atmospheric  emission  or  background  radiation. 

Error  is  also  introduced  due  to  the  neglect  of  temperature  dependence  in  the 
model.  The  spectral  curves  of  Figs.  16-25  were  constructed  for  a  temperature 
of  27 3K  and  atmospheric  temperatures  either  cooler  or  warmer  provide  some 
variability.  This  effect  is  not  too  large  for  normal  atmospheric  temperatures 
except  in  the  extreme  wings  of  bands  and  in  bands  arising  from  an  excited  lower 
energy  state  where  the  population  of  the  energy  Btates  is  quite  temperature 
dependent..  This  is  usually  only  true  for  weak  bands  or  in  weak  portions  of  stronger 
absorption  bands  and  so  is  not  too  critical  for  our  purposes.  In  general,  cooler 


temperatures  will  mean  somewhat  higher  transmittance  values  and  higher  tempera¬ 
tures  somewhat  lower  transmittance  values. 

It  cannot  be  emphasized  enough  that  the  results  obtained  from  these  low  spectral 
resolution  curves  should  never  be  used  in  the  determination  of  atmospheric  attenua¬ 
tion  for  laser  propagation.  Detailed,  very  high  spectral  resolution  measurements 
or  calculations  such  as  are  presented  m  Table  5  and  elsewhere  (IVicClatchey,  1971 
and  McClatchey  and  Selby,  1972)  are  required  for  this  purpose. 

6.  SCATTERED  SOLAR  RADIATION  AND  INFRARED  EMISSION 

A. I  Transmitted  and  Reflected  Sky  Radiance  in  Clear  and  Cloudy  Sky 

While  the  energy  taken  out  of  a  beam  of  radiation  by  absorption  contributes  to 
heating  of  the  air,  the  energy  scattered  by  molecules,  aerosols,  or  cloud  droplets 
will  be  redistributed  in  the  atmosphere.  We  will  only  consider  the  problem  of  sun¬ 
light  incident  on  a  plane  parallel  atmosphere.  As  long  as  only  single  scattering  is 
important  and  a  homogeneous  atmosphere  is  assumed  (aerosol/air  mixing  ratio 
constant  with  height),  the  angular  distribution  of  skylight  at  a  specific  wavelength 
can  be  obtained  from 


B(M,  =  r 


i  m'  \ 

/  -tM*  -tM) 

\M  -  M'J 
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1  =  extraterrestrial  solar  irradiance  (W  cm 

o 

-2  -1 

B  =  sky  radiance  (W  cm  sterad  ) 

M  =  secant  of  the  solar  zenith  angle  6 

M1  =  secant  of  the  line  of  sight  angle  e' 

t  -  a  (0  +  k  )  -t-  m  (a  +  k_)  -  total  optical  thickness  per  unit  air  mass 

V  d  <1  V  ill  fit 

(scattering  and  absorption  by  molecules  and  aerosola-see  Section  5.3 

for  ay  and  mv;  <ra,  k^,  and  crm  can  be  derived  from  Figures  21  and  22 

by  taking  the  logarithm  of  the  transmittance).  With  regard  to  molecular 

absorption  (k  ),  Eq.  24  is  only  valid  for  monochromatic  radiation 

(sec  Table  5  for  k  values) 
m 

t  *  a  c  +  m  a  =t  +  t 
s  v  a  v  m  s,a  s,m 

FW>)  =  angular  scattering  function  (molecular  +  aerosol)  per  unit  air  mass 
=  ($)  t  +  t  (see  Figure  26  for  <$?(<#>)  and 

ci  o  it*  Sjtri  3  jit 
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The  scatter  ing  angle  </b  is  defined  by  cos  <f>  -  cos  0  cos  O'  -  sin  0  sinG1  cos  (A'l/) 
where  A'}'  is  the  angular  azimuth  difference  between  the  sun  direction  and  the  line 
of  sight. 

However,  if  tM  or  tM '  becomes  larger  than  0.10,  higher  order  scattering  can 
no  longer  be  neglected;  this  means  that  illumination  of  the  scattering  volume  from 
the  sky  and  reflecting  earth's  surface  becomes  increasingly  important.  Higher 
order  scattering  becomes  dominant  for  tM  or  tM1  >0.5,  particularly  for  high 
ground  albedo.  This  is  particularly  true  for  the  radiation  leaving  the  top 
of  the  atmosphere  (radiation  to  space).  The  general  radiative  transfer 
equation  has  been  solved  analytically  only  for  molecular  scattering, 
and  detailed  tables  of  sky  radiance,  polarization  (including  the  Stokes  parameters) 
have  been  prepared  by  Coulson  et  al<1960).  Skylight  distributions  for  hazy  atmos¬ 
pheres  have  been  computed  by  de  Bary  et  al  (1965)  using  these  tables  and  Eq.  (24). 
However,  computers  have  made  possible  a  statistical  treatment  of  the  scattering 
problem  in  any  realistic  hazy  or  cloudy  atmosphere.  By  this  Monte  Carlo  method, 
several  investigators  have  computed  the  distribution  of  the  radiance  of  the  hazy 
atmosphere,  both  as  seen  from  the  ground  (Figure  27a)  and  from  space  (Figure  27b). 
The  assumed  aerosol  size  distribution  and  vertical  distribution  is  essentially  the 
same  as  in  Table  3.  Dotted  curves  of  the  radiance  (per  unit  solid  angle  and  unit 
incident  solar  flux)  are  for  the  sun  in  the  zenith,  and  solid  lines  are  for  the  sun  at 
6  =  86.  3°. 

In  Figure  27a,  the  downward  or  "transmitted"  radiance  generally  increases 
with  decreasing  wavelength  when  the  sun  is  in  the  zenith.  Forward  scattering 
causes  the  radiance  to  peak  near  the  sun.  Near  the  horizon,  the  radiance  increases 
again,  except  at  short  wavelengths,  and  the  albedo  influence  is  large.  At  low  solar 
elevations,  downward  radiance  is  generally  much  smaller,  especially  near  the 
zenith  at  long  wavelengths,  and  the  albedo  influence  is  significant  only  at  short  wave¬ 
lengths.  Since  the  sky  radiance  values  have  been  averaged  over  intervals  of  0.  1 
in  cob  8  (that  is,  over  angular  intervals  ranging  from  5°  to  25°),  the  radiance  near 
the  sun  actually  should  be  much  higher  than  shown. 

In  Figure  27(b),  upward  or  reflected  radiance  as  seen  from  outside  of  the 
atmosphere  is  strongly  dependent  on  the  albedo  of  the  ground.  The  contribution  of 
reflected  ground  radiance  becomes  dominant,  especially  in  the  near  1R  and  in  the 
nadir  direction  where  the  atmospheric  backseat^ erect  flux  becomes  umall.  This 
effect  will  be  less  pronounced  in  more  hazy  atmospheres. 

The  distribution  of  sky  radiance  in  the  ultraviolet  is  dominated  by  multiple 
scattering  and  below  0.  35  pm  by  ozone  absorption.  At  high  colar  elevations,  down¬ 
ward  as  well  as  upward  radiances  at  0.  30  pm  are  only  about  1  percent  of  their  res¬ 
pective  values  at  0.  40  pm,  but  independent  of  albedo- 

In  the  case  of  the  radiance  from  clouds,  we  present  Figures  28a  and  28b  for 
a  rather  dense  nimbostratus  (Kattawar  and  Plass,  1968).  Aerosol  and  molecular 
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scattering  have  been  neglected  in  this  model  because  they  would  have  little  effect. 
The  strong  downward  radiance  (Figure  2  8a)  for  high  sun  and  low  albedo  is  largest 
near  the  sun,  but  with  high  albedo  is  nearly  independent  of  the  angle  of  observation. 
At  low  sun,  the  much  smaller  cloud  radiance  is  largest  in  the  zenith.  The  radiance 
of  the  cloud  top  (Figure  28b)  again  is  very  large  at  high  sun.  At  low  sun  there  is 
little  discrimination  between  a  cloudy  and  a  clear  sky  (see  Figure  27,  \  =  0.  7  pm). 

An  indication  of  the  variability  of  the  reflectance  or  spectral  albedo  of  some 
natural  objects  may  be  obtained  from  Figure  29.  The  data  presented  are  typical 
with  respect  to  the  waveiengih  dependence,  but  absolute  reflectance  may  vary  con¬ 
siderably  for  individual  objects,  and  generally  is  higher  for  radiation  near  grazing 
incidence. 

6.2  Infrared  Emission  of  the  Atmosphere 

In  the  infrared  from  5  to  25  pm  and  beyond,  the  atmospheric  emission  in  a 
collimated  beam  can  be  determined  by  numerical  integration  of  the  equation  of 
radiative  transfer  to  give  the  upward  radiance,  IT  (v): 


11  (v) 


-X 


B  (v,  T)  dT  +  B  (v,  T  )  r 

O  O 


where  the  integral  represents  the  atmospheric  contribution  and  the  second  term  is 
the  contribution  from  the  ground  (or  lower  boundary)  modified  by  the  total  atmos¬ 
pheric  transmittance  t  between  the  bqundary  and  the  observer;  r  is  the  transmit- 
6 

tance  which  ranges  between  1.  0  at  the  level  of  observation  and  at  the  lower 
boundary  (for  example,  the  earth's  surface  or  cloud  top);  B  (v,  T)  is  the  Planck 
(black  body)  function  corresponding  to  frequency  v  and  the  temperature  T  of  an 
atmospheric  layer,  and  Tg  is  the  ground  temperature. 

To  obtain  the  downward  emission  from  a  clear  sky,  the  second  term  in  Eq.  (25) 
should  be  omitted  and  Ih  then  the  transmittance  of  the  entire  atmosphere.  In  a 
cloudy  sky  the  subscript,  g,  refers  to  a  cloud  base. 

The  effects  of  aerosols  and  of  molecular  scattering  nave  been  neglected  and 
the  molecular  absorption  charts  described  in  Section  5. 2  have  been  applied  to  the 

water  vapor,  ozone,  and  uniformly  mixed  gases  corresponding  to  the  Midlatitude 
Winter  model.  Equation  (25)  was  then  applied  using  the  derived  transmittance 
values  and  performing  the  numerical  integration  according  to  a  numerical  analogue 
of  Eq.  (25). 

Ihe  results  are  presented  in  Figure  30  for  both  downward  and  upward  beam 
radiation.  The  downward  radiation  case  assumes  the  observer  to  be  at  the  ground 
looking  up  at  the  zenith.  The  upward  case  assumes  an  observation  point  at  the  top 
of  the  atmosphere  viewing  the  nadir  (that  is,  looking  down). 
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It  should  be  noted  that  in  the  upward  case  the  major  absorption  features  (6.  3jum 
H2O,  9.6  pm  O3,  and  15  Aim  CO2)  appear  as  relative  minima,  whereas  in  the  down¬ 
ward  case  these  same  regions  appear  as  relative  maxima.  The  surface  tempera¬ 
ture,  Tg,  was  taken  to  be  the  same  as  the  temperature  of  the  lowest  atmospheric 
level  in  the  model  (272K).  The  upper  curve  in  Figure  30  corresponds  to  a  black- 
body  at  this  boundary  temperature. 

For  wavelengths  shorter  than  5  Aim,  the  problem  of  upward  and  downward 
radiance  is  complicated  by  molecular  and  aerosol  scattering  and  is,  therefore,  a 
function  of  the  solar  zenith  angle.  Computations  of  radiation  between  3  and  5 
micrometers  is  particularly  complicated.  Some  examples  of  measured  radiation 
in  the  1  to  5  pm_  region  are  included  in  Figure  31. 


7.  REKRU  TIVK  K1TTXTS  IN  THE  ATMOSPHERE 


Previously  (Sections  3  and  5)  it  was  indicated  that  slant  path  length  could  be 
obtained  from  vertical  path  length  through  multiplication  by  sec  0  for  0  <  80°.  For 
0  2  80°,  the  vertical  path  length  (or  other  concentration  parameter),  AL,  may  be 
converted  to  slant  path  length  by 


AL  (slant  path)  =AL  (Vertical)  x  M(0) 


(26.) 
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where  M(0)  is  the  relative  optical  path,  defined  as  the  ratio  of  slant  path  length  to 
vertical  path  length  between  two  levels.  The  term  M(0)  always  equals  sec  0  for  a 
short  path  near  the  observer,  and  for  any  path  if  0  <  80°.  For  a  long  slant  path, 
however,  M(0)  becomes  smaller  than  sec  0  because  of  the  curvature  of  the  atmos¬ 
phere,  which  is  partly  compensated  by  refraction  of  tire  beam.  This  implies  that 
the  vertical  distribution  of  a  scattering  or  absorbing  atmospheric  constituent  also 
influences  the  optical  path,  and  the  slant  path  length  may  have  to  be  derived  from 
Eq.  (26)  if  nigh  accuracy  is  warranted. 

If  curvature  and  refraction  are  considered,  the  genera)  equation  for  the  optical 
path  between  altitudes  zq  and  z  ^  of  a  constituent  of  density  p  at  z  and  scale  height 


11  =  (Zj  -  zQ)  /logc  [  p(z0)/p(z1) ]  may  be  written 

M(0)  =[{W0>  '  P(zi>)H]  1  fz  [X  "(7+r)  ("Sfe  Sin  6)2] 


p(z)  dz 


(27) 


where  n  is  the  refractive  index  of  air  and  the  (o)  subscript  refers  to  the  lower 
leference  level,  and  r  is  the  radius  of  the  earth. 
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For  the  molecular  atmosphere  and  uniformly  mixed  atmospheric  constituents, 
the  air  mass  relative  to  the  vertical  air  mass  of  the  entire  atmosphere  measured 
from  the  ground,  M(0),  as  computed  by  Kasten  ( 1967 )  is  presented  in  Figure  32. 
These  data  are  based  on  the  ARDC  Model  Atmosphere  of  Minzner  et  al  (1959)  and 
a  refractive  index  of  air  corresponding  to  a  wavelength  of  0.7  pm.  Earlier  values 
of  M(9)  of  which  the  data  by  Bemporad  (1907)  are  most  often  used,  are  in  fair 
agreement.  Kasten  (1967)  has  shown  that  the  average  influence  of  latitude  and 
season  on  M(0)  is  less  than  2  percent  but.  abnormal  tropospheric  density  profiles 
may  have  larger  effects.  Figure  32  compares  M(0)  values  for  dry  air  with  Mw(0) 
values  for  an  atmosphere  containing  water  vapor  with  a  distribution  according  to 
Schnaidt  (1938).  Values  of  sec  0  are  included  for  comparison 

The  angular  deviation  between  the  actual  and  apparent  solar  zenith  angles,  AO, 
due  to  refractive  effects  is  presented  in  Table  C. 


Table  6.  Difference  between  Actual  and  Apparent  Zenith  Angle  as  a  Function 
of  Apparent  Zenith  Angle 


Vp(DEG) 

75 

80 

83 

85 

86 

87 

A0  (DEO) 

0.  06  i 

0.  080 

0.  124 

0.  167 

0.  199 

0.248 

%p  <DEG> 

88 

88.  5 

89 

89.  5 

90 

A0  (DEG) 

0.  311 

0.  358 

0.417 

0.497 

0.  570 

Table  7.  Relative  Optical  Pathol'  a  Thin,  High  Altitude  Layer  as  a 
Function  of  Zenith  Angle  0 


Height  of  Layer 

(km)  Above 

Sui'faCC 

0 

10 

20 

30 

40 

50 

60 

70 

2.  89 

2.  86 

2.  83 

2.  60 

2.  77 

2.  74 

75 

3.78 

3.71 

3.  63 

3.  57 

3.  50 

3.44 

80 

5.49 

5.26 

5.05 

4.  87 

4.70 

4.  56 

85 

9.  66 

8.  51 

7.70 

7.  OB 

6.  60 

6.  20 

87 

13.07 

10  56 

9.  1) 

8.  13 

7.41 

6.  86 

88 

15.  17 

11.  58 

9.74 

8.  57 

7.74 

7.  1? 

89 

17.  08 

12.  37 

10.  19 

8.  87 

7.96 

7.28 

89.  5 

17.  66 

12. 58 

10.  30 

8.  94 

8.  01 

7.  33 

1  90 

17.88 

12.  66 

10.  35 

8.97 

8.03 

7. 34 
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For  thin  atmospheric  layers  (10  km  or  less)  at  high  altitudes  (where  refractive 
effects  can  be  neglected).  Table  7  provides  the  geometric  portion  (that  is,  the  ef¬ 
fect  is  due  to  curvature  only)  of  the  slant  path  correction.  For  example,  the  geo¬ 
metric  path  through  a  thin  layer  at  30  km  viewed  at  a  zenith  angle  of  90°  is  10.  36 
times  the  geometric  path  in  the  vertical.  However,  the  optical  path  through  a 
uniformly  mixed  atmospheric  constituent  at  the  same  viewing  angle  would  be  36.  3 
(see  Figure  32).  This  situation  may  occur  in  the  treatment  of  transmittance  through 
the  ozone  layer  or  stratospheric  dust  layers. 


8.  KHTXT  01'  CI.OIDS  \ND  HX,  ON  ATMOStMlKHIO  RADIATION 

As  in  the  case  of  aerosol  extinction  (Section  3.  2)  when  a  beam  of  radiation  is 
intercepted  by  a  cloud  or  fog  layer,  some  of  the  incident  radiation  may  be  absorbed 
and  some  may  be  scattered  out  of  this  beam  with  both  effects  thereby  causing  a  net 
attenuation  of  the  original  radiation.  The  magnitude  of  the  attenuation  depends  on 
the  number,  size  distribution  and  refractive  index  of  the  suspended  droplets  and 
on  the  thickness  of  the  layer.  Since  the  radiation  may  be  scattered  in  all  direc¬ 
tions,  including  those  close  to  that  of  the  original  beam,  the  resulting  attenuation 
may  also  depend  on  the  acceptance  angle  of  the  receiver. 

Natural  fogs  and  low  level  clouds  are  composed  of  spherical  water  droplets, 
the  refractive  properties  of  which  have  been  fairly  well  documented  in  the  spectral 
region  of  interest  here  (Kislovskii,  1959;  and  Irvine  and  Pollack,  1968).  However, 
in  the  case  of  cirrus  clouds  the  particles  consist  of  ice  crystals  of  various  shapes 
(hexagonal  or  cylindrical  crystals)  which,  when  mutually  oriented,  can  give  rise 
to  anomalous  scattering  and  polarization  phenomena  when  radiation  is  propagated 
through  them. 

In  this  section  we  will  restrict  outseives  to  the  propagation  of  a  collimated 
beam  of  radiation.  Generally,  fogs  and  clouds  show  only  a  small  increase  in  trans¬ 
mittance  with  wavelengths  from  the  visible  to  about  6  pm.  Beyond  6  pm,  a  greater 
increase  is  observed.  However,  certain  infrequently  occurring  types  of  fog 
characterized  by  a  sharply  peaked  distribution  of  small  droplets  (this  phenomenon 
is  also  found  near  the  base  of  stratus  and  cumulus  clouds),  show  a  significant  in¬ 
crease  in  infrared  transmittance. 

These  general  effects  can  be  seen  from  Figure  33,  which  shews  the  wavelength 
variation  of  the  extinction  coefficient  for  a  drop  size  distribution  representative  of 
fair  weather  cumulus  cloud  after  Deirmendjian  (1964  and  1969).  The  model  drop 
size  distribution  (shown  in  Figure  34),  which  is  given  by 

n(r)  =  2.  373  r®  e”*"  ®r  cm  ®  pm 


(28) 
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has  a  maximum  corresponding  to  a  radius  of  4  pm  and  is  normalized  to  100  drops 
cm"3 

Figure  33  was  calculated  and  based  on  single  scattering  theory,  using  the 
complex  refractive  index  values  for  water  given  by  Kislovskii  (1959)  and  Eq.  (29). 

It  will  be  seen  from  Figure  33  that  the  attenuation  due  to  this  fairly  thin  model 
cloud  is  rather  severe;  for  example,  for  a  100m  path  length  the  transmittance  is 
approximately  0.  3  and  0.  15  at  10  pm  and  1  pm  respectively. 

Some  mention  should  be  made,  at  this  point,  of  the  validity  of  using  single 
scattering  theory  to  account  for  the  attenuation  of  radiation  in  cloud. 

It  has  generally  been  assumed  that  calculations  of  attenuation  based  on  single 
scattering  theory  are  valid  for  an  optical  thickness  of  less  than  0.03  (van  do  Hulsts 
criterion)  and  that  for  values  of  optical  thickness  in  excess  of  unity,  secondary  and 
multiple  scattering  effects  become  important.  However,  this  applies  to  the  propa¬ 
gation  of  diffuse  radiation.  For  a  narrow  beam  of  collimated  radiation  (or  laser 
beam  propagation)  it  has  been  found  (Zuev  et  al,  1967)  that  single  scattering  theory 
can  account  for  the  observed  attenuation  for  values  of  optical  thickness  of  up  to  25. 
In  this  context  optical  thickness  is  defined  as  the  product  of  the  extinction  coeffi¬ 
cient  and  path  length. 


Figure  1.  Normalized  Particle  Size  Distribution 
for  Aerosol  Models 


Figure  2.  Equivalent  Sea  Level  Path  Length  of  Water  Vapor  as  a 
Function  of  Altitude  for  Horizontal  Atmospheric  Paths 
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Figure  4.  Equivalent  Sea  Level  Path  Length  of  Uniformly  Mixed  Gases 
as  a  Function  of  Altitude  for  Horizontal  Atmospheric  Paths 
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Figure  5.  Equivalent  Sea  Level  Path  Length  of  Uniformly  Mixed  Gases 
as  a  Function  of  Altitude  for  Vertical  Atmospheric  Paths 
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Figure  12.  Equivalent  Sea  Level  Path  Length  for  Molecular  Scattering 
as  a  Function  of  Altitude  for  Horizontal  Atmospheric  Paths 


Figure  16a.  Prediction  Chart  for  Water  Vapor  Transmittance  (0.  6-4.  Oum) 
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Prediction  Chart  for  Ozone  Transmittance  <3-6.4  fim) 
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Fijfure  13b.  Prediction  Chart  for 
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Prediction  Chart  for  Nitrogen 
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Prediction  Chart  for  Transmittance  Due  to  Water  Vapor  Continuum 
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for  Transmittance  Due  to  Molecular  Scatte: 


Chart  for  Aerosol  Transmittance  (Scattering  and  Absorption) 


Ozone  (0.  25-0.  75  ym)  as  a  Function  of  Altitude 
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ure  24.  Equivalent  Sea  Level  Path  Length  Tor  Ozone  (0.  2 5-0.  75  »im)ass  Function  of  Altitude 
Vertical  Atmospheric  Paths 
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Figure  31.  Spectral  KaHiance  of  3  Clear  Xenith  Sky  Showing  the  I'ffect  of  the 
Position  of  the  Sun  on  the  (Initiation  Scattered  Below  3  anil  the  Kffect  of 
Thermal  emission  at  linger  Wavelengths  Measured  at  Cocoa  Beach, 
Florida  (see  Bell  et  al,  1  Of >0 ) 
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Figure  32.  Comparison  of  the  r.clativr  Air  Mass  [M«?)]Oomputcd  by  Hasten  ll%7) 
as  a  Function  of  Apparent  Zenith  Angie.  6,  llasid  on  the  ARDC  Model  Atmosphere 
(Minznor  ot  al.  1999)  With  a  Similar  (Quantity  Rased  on  a  vertical  Water 

Vapor  Distribution  According  to  Schnaidt  (1938)  and  for  sec  6  (Note  that  for  ©<00°, 
m (0)  "  Mw«?)  -  sic©.) 
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Function  of  Wavelength  Due  to  Fair-weather  Cumulus  Cloud 


Distribution  for  Fair-weather  Curr.ulus  Cloud  Model 
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Appendix  A 

Units  and  Conversion  Factors 

1.  INTRODUCTION 

Throughout  this  report  we  have  used  ihe  term  "frequency"  rather  loosely  to 
signify  spectral  variation.  In  general,  frequency  is  assumed  to  be  synonymous 
with  wavenumber,  v  (cm'*),  which  is  given  by  v  =  104/A  where  the  wavelength  X 
is  given  in  units  of  p rn  (that  is,  10“^  meter). 

The  subject  of  units  used  in  atmospheric  physics  is  rather  confusing  to  the 
beginner.  This  is  particularly  true  where  units  of  absorber  concentration  or 
amount  are  concerned.  In  this  Appendix  we  will  show  the  reader  how  to  derive  the 
absorber  concentrations  for  use  with  the  prediction  charts  given  in  Section  5,  on 
the  assumption  that  one  has  the  relevant  meteorological  data  available.  We  will 
then  discuss  the  units  of  absorber  concentration  and  indicate  how  to  make  the  ap¬ 
propriate  conversions. 

2.  ABSORBER  CONCENTRATION 

For  most  gases  absorber  concentration  denoted  here  by  (S L  is  generally  mea¬ 
sured  by  the  equivalent  length  of  pure  absorber  in  the  path  under  the  existing  condi¬ 
tions  of  pressure,  P,  and  temperature,  T.  The  unit,  centimeter  atmospheres, 

(cm -atm)  is  given  by  &L  s  cRP  where  c  is  the  fractional  concentration  of  the  ab¬ 
sorber  (Bee  Table  2),  R  is  the  path  length  (cm)  and  the  P  is  the  pressure  in  atmos¬ 
pheres.  More  often  AL,Q  is  quoted  as  the  absorber  concentration  reduced  to  GTP 
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conditions,  which  is  a  more  meaningful  quantity  where  vertical  or  slant  atmospheric 
paths  are  concerned.  Thus,  for  a  horizontal  path  at  altitude,  z,  we  have 


AL 


o 


:cR 

O 


(tF>)  <<™-*tm>STP 


(A.  1) 


whe-  *  273K,  and  PQ  =  1  atm.  If  c  is  given  in  parts  per  million  (ppm),  (see 
Tab*  2),  and  R  is  given  in  km.  the  following  expression  can  be  used. 


ALo  (cm-atm)g^,p  =  10 


-1 


•■*(*) 


<A2) 


For  a  uniformly  mixed  gas  of  fractional  concentration  by  volume,  c,  the  number  of 
centimeter  atmospheres  at  STP  in  a  slant  path  from  any  altitude  z  to  the  limit  of 
the  atmosphere  is  given  by 


ALo  (0,  z 


»  ), 


STP 


c  H  (z)  sec  6 


p(z)  (  To  \ 

Pn  \W)/ 


*  c  H(zl  sec  0  p(z)  /  p 


(A3) 


where  H(z)  is  called  the  scale  height  of  the  atmosphere  above  altitude  z,  6  is  the 
zenith  angle  and  p(z)  and  pq  are  the  atmospheric  density  at  altitude  z  and  STP. 
respectively.  For  6  >  70°.  see  Section  7.  From  the  hydrostatic  equation,  H(z) 
can  be  defined  as  follows: 


P(z)  =  g  p  (z)  H(z) 


(A4) 


where  P(z)  is  the  pressure  at  altitude  z,  and  g  is  the  acceleration  due  to  gravity 
(981  cm  sec~2).  Thus,  for  a  vertical  column  of  air  at  STP.  H(o)  *  7.  99  x  10°  cm 
(that  is.  approximately  3  kin).  If  the  scale  height  in  cm  is  required  at  another 
altitude,  z.  it  can  be  determined  from 


H(z)  *  1.02  x  106  x  iF(z)/p  <z)| 


(A  5) 


where  P(z)  is  given  in  millibars  and  p(z)  iB  given  in  gm/m^  (see  Table  1).  Note 
from  Eq.  (A 5)  that  the  quantity  H(z)  p(z)  is  directly  proportional  to  P(z)  and  is 
independent  of  temperature.  Thus,  inserting  H(z)  p(z)  into  Eq.  (A3),  we  have 


ALq  (6,  z  —  <«  )STp  *  c  x  1.  02  x  10S  P(z>  sec  6  /pQ 


(A6) 


{ 


S  and  for  the  amount  of  absorber  between  levels  Zj  and  (see  Table  2): 

i  ALo  (e.Zj- Z2)STP  =  e  x  1.02  x  106  [  PCz.j)  -  P(z2))  sec  e/pQ.  <A7) 

For  condensible  vapors  (for  example,  water  vapor)  it  is  more  convenient  to 
measure  the  total  length  of  liquid  (in  ems)  which  may  be  precipitated  out  of  the  path 
(per  unit  area).  The  number  of  precipitable  centimeters  (abbreviated  pr  cm)  of 
water  vapor  in  a  given  path  of  length  R(cm)  is  given  by .  w  -  p(T)  WR  where  p(T)  is 
the  saturation  vapor  density  of  water  at  temperature  T  (gm  cm-3)  and  W  is  the 
fractional  relative  humidity.  It  will  be  noted  that  the  units  of  pr  cm  and  gm  cm~2 
are  equivalent. 

Figure  A1  will  enable  the  reader  to  calculate  the  number  of  pr  ems  of  water 
vapor  per  kin  path  length,  for  a  given  temperature  and  relative  humidity.  For 
temperatures  below  -40°C  the  reader  should  use  density  values  of  water  vapor 
measured  over  ice  for  determining  water  vapor  concentrations  for  atmospheric 
paths  (see  Smithsonian  Meteorological  Tables.  1963). 

In  Table  1,  water  vapor  densities  are  given  in  units  gm  m  '  as  a  function  of 
altitude  for  the  5  model  atmospheres.  It  may  be  useful  to  remember  that  fo’-  water 
vapor . 

1  gm  m'3  =  0.  1  gm  cm'2/km  (or  0.  1  pr  cm/km). 


3.  (  NIT  (  OWKRSION  FACTORS 

According  to  Avogadros  hypothesis  the  molecular  weight,  M.  in  gms.  ,  or  any 
gas  occupies  22. 4  liters  at  STP.  Hence  one  cubic  centimeter  of  gas  at  STP  weighs 
(M/2.24.  104)  gms.  Since  1  (cm-atm)spp  is  equivalent  to  a  length  of  1  cm  of  gas 
at  STP  per  cm3,  we  can  write  the  following: 


1  cm  -atm  of  gas 


M 


v  in4 


gm  cm 


(A8) 


where  M  is  the  molecular  weight  of  the  gas  (see  Table  2). 

Also  since  1  cm^  of  gas  at  STP  contains  1/2.  24x  104  tnoles  and  one  mole  of  gas 
contains  the  same  number  (ft  )  of  molecules  (where  &  is  Avogadros  Number)  we 
can  write  that  1  cm3  of  gas  at  STP  contains  6.  06 x  10^3/2, 24  x  104  *  2.  69x  10 ^ 
molecules  or 


1  (cm-atm)gTp  =  2.69x  10  ly  molecules/cm2. 


(A9) 


<**»**»* 


!W 


Note  that  Kq.  (AD)  is  valid  independent  of  the  gas.  Since  for  water  vapor  one  can 
write  I  g;n  cm'2  in  terms  of  (cm -atmj^.pp,  using  Kq.  (A8)  it  will  also  be  seen  that 
for  water  vapor: 

1  pr  cm  (H.jO)  =  3.34x  10^  molccules/cm'^ .  (A10) 

Thus  the  unit,  moleculc/cm*2 ,  is  independent  of  the  nature  of  the  absorbing  gas  and 
basic  to  all  gases.  It  is  hoped  that  this  unit  will  become  more  widely  used  and 
thereby  help  to  unify  absorber  concentration  units  for  all  the  atmospheric  cor. 
tuents. 
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Appendix  B 

Equivalent  Sco  Level  Path  Charts  for  U.S.  Standard  Atmosphere,  1962 


In  both  the  original  version  of  this  report  (McClatchcy,  et  al,  1970)  and  in 
the  revised  version  (McCiatchey,  et  al,  1971),  it  was  decided  to  provide  explicit 
computational  capability  for  five  atmospheric  molecular  models  and  two  aerosol 
models.  The  inclusion  of  water  vapor  models  in  the  U.S.  Standard  Atmosphere 
Supplements,  1906  as  provided  in  the  2  FCKL  Handbook  of  Geophysics  and  Space 
Environments  (Valley,  1965),  provided  an  appropriate  set  of  models  for  this  pur¬ 
pose.  The  U.S.  Standard  Atmosphere,  1962  does  not  provide  an  accompanying 
water  vapor  distribution  and  so  this  model  was  initially  omitted. 

Within  the  meteorology  and  atmospheric  physics  community,  there  is  a  need 
to  intercompare  calculations  for  a  single  appropriate  atmospheric  model.  It  is 
felt  that  the  best  model  for  this  purpose  Is  the  U.S.  Standard  Atmosphere,  1962 
with  an  appropriate  water  vapor  and  ozone  distribution.  The  U.S.  Standard 
Atmosphere  Supplements,  1966  provides  an  appropriate  ozone  distribution  utilizing 
the  experimental  work  of  Bering  and  Borden,  1964  and  calculations.  The  model 
specified  there  is  given  below  in  Table  Bl.  An  appropriate  water  vapor  distribution 
(Sissenwine,  1961)  has  also  been  included  in  Table  Bl  together  with  the  tempera¬ 
ture,  pressure  and  density  results  from  the  U.S.  Standard  Atmosphere,  196!. 

Using  the  model  atmosphere  parameters  from  Table  Bl  and  the  equivalent 
sea  level  path  relations  from  Section  5.  3,  Figures  Bl  -  B14  have  been  constructed. 
These  figures  should  be  used  to  determine  the  equivalent  sea  level  path  lengths 
for  the  Standard  Atmosphere  and  the  results  used  with  Figures  16a  through 
Figure  25  as  described  in  Section  5.4. 


Table  Bl.  Model  Atmosphere  Used  as  a  Basis  for  the  Computation 
of  Atmospheric  Optical  Properties 


U.  S.  STANDARD  ATMOSPHERE,  1962 


Ht. 

(km) 

Pressure 

(mb) 

Temp. 

<°K) 

Density 
(g  /  m3) 

Water  Vapor 
(g/m3) 

Ozone 

(g/m3) 

0 

1.  013E+03 

288.  1 

1.  225E+03 

5.  9E+00 

5.  4E-05 

1 

8.  986EH-02 

281.  6 

1.  111E+03 

4.  2E+00 

5.  4E-05 

2 

7.  950E+02 

275.  1 

1.  007E+03 

2.  9E+00 

5.  4E-05 

3 

7,  012E+02 

268.  7 

9.  093E402 

1.  8E+00 

5.  0E-05 

4 

6.  1C6E+02 

262.  2 

8.  193E+02 

1.  1E+00 

4.  6E-05 

5 

5.  405E+02 

255.  7 

7.  364E+02 

6.  4E-01 

4.  5E-05 

6 

4. 722E+02 

249.  2 

G.  601E+02 

3.  8E-01 

4.  5E-05 

7 

4.  111E+02 

242.  7 

5.  900E+02 

2.  IE-01 

4.  8E-05 

3 

3, 565E+02 

236.  2 

5.  258E+02 

1.  2E-P1 

5.  2E-05 

9 

3.  080E102 

229.  7 

4.  671E+02 

4.  6E-02 

7.  IE-05 

10 

2.  650E+02 

223.  2 

4.  135E402 

1.  8E-02 

9.  9K-05 

11 

2. 270E+02 

216.  8 

3.  648E+02 

8.  2E-03 

1.  3E-04 

12 

1.  940E+02 

2 16.  6 

3.  119E+02 

3.  7E-03 

1.  6E-04 

13 

1, 658E+02 

216.  6 

2.  666E+02 

1.  8E-03 

1.7E-04 

14 

1.  417E+02 

216.  6 

2.  279EH  02 

8.  4K-04 

1.  9E-04 

15 

1.  211E+02 

216.  6 

1.  948E+02 

7.  2E-04 

2.  IF.- 04 

16 

1.  O35E+02 

216.  6 

1.  665E+02 

6.  IE- 04 

2.  3E-04 

17 

8. 850E+01 

216.  6 

1.  423E+02 

5.  2E-04 

2.  BE-  04 

18 

7. 565E+01 

216.  6 

1.  216E+02 

4.  4E-04 

3.  2E-04 

19 

6.  467E+01 

216.  6 

1.  040E+02 

4.  4E-04 

3.  5E-04 

20 

5.  529E+01 

216.  6 

8,  891E+01 

4.  4E-04 

3.  8E-04 

21 

4.  729E+01 

217.  6 

7.  572E+01 

4.  8E-04 

3.  8E-04 

22 

4,  047E+01 

218.  6 

6.  45iE+Gi 

5.  2E-04 

3.  9E-04 

23 

3.  467E+01 

219.  6 

5.  500E+01 

5.  7E-04 

3.  8E-04 

24 

2.  972E+01 

220.  6 

4.  694E+01 

6. IE-04 

3.  6E-04 

25 

2.  549E+01 

221.6 

4.  008E+01 

6.  6E-04 

3,  4E-04 

30 

1.  197E+01 

226.  5 

1. 841E+01 

3.  8E-04 

2.  0E-04 

35 

5.  746E+00 

236.  5 

8.  463E+00 

1.  SE-04 

1.  IE-04 

40 

2.  871E+00 

250.  4 

3.  906E+00 

6.  7E-C5 

4.9E-05 

45 

l. 491E+00 

264.  2 

1.  966E+00 

3,  2E-05 

1.7E-05 

50 

7. 978E-01 

270.6 

1.  027E+00 

1.  2E-05 

4.  0E-G6 

70 

5. 520E-02 

219.  7 

8.  754E-02 

1.  5E-07 

8.  6E-08 

100 

3. 008E-04 

210.  0 

4.  989E-04 

1.  0E-09 

4.  3E-11 

ALTITUDE  (km) 


Figure  Bl.  Equivalent  Sea  Level  Path  Length  of  Water  Vapor  as  a  Function  of 
Altitude  for  Horizontal  Atmospheric  Paths 
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Figure  B2.  Equivalent  Sea  Level  Path  Length  of  Water  Vapor  as  a  Function  of 
Altitude  for  Vertical  Atmospheric  Paths 
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Figure  B3.  Equivalent  Sea  Level  Path  Length  of  Uniformly  Mixed  Gases  as  a 
Function  of  Altitude  for  Horizontal  Atmospheric  Paths 
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Figure  BIO.  Equivalent  Sea  Level  Path  Length  for  Water  Vapor  Continuum  as 
a  Function  of  Altitude  for  Vertical  Atmospheric  Paths 
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Figure  Bll.  Equivalent  Sea  Level  Path  Length  for  Molecular  Scattering  as  a 
Function  of  Altitude  for  Horizontal  Atmospheric  Paths 
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Figure  B12.  Equivalent  Sea  Level  Path  Length  for  Molecular  Scattering  as  a 
Function  of  Altitude  for  Vertical  Atmospheric  Paths 
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